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We present final results on the photon electroproduction {ep —^ ep 7 ) cross section in the deeply 
virtual Compton scattering (DVCS) regime and the valence quark region from Jefferson Lab ex¬ 
periment EOO-110. Results from an analysis of a subset of these data were published before, but 
the analysis has been improved which is described here at length, together with details on the 
experimental setup. Furthermore, additional data have been analyzed resulting in photon electro¬ 
production cross sections at new kinematic settings, for a total of 588 experimental bins. Results of 
the Q^- and xs-dependences of both the helicity-dependent and helicity-independent cross sections 
are discussed. The Q^-dependence illustrates the dominance of the twist-2 handbag amplitude in 
the kinematics of the experiment, as previously noted. Thanks to the excellent accuracy of this 
high luminosity experiment, it becomes clear that the unpolarized cross section shows a significant 
deviation from the Bethe-Heitler process in our kinematics, compatible with a large contribution 
from the leading twist-2 DVCS^ term to the photon electroproduction cross section. The necessity 
to include higher-twist corrections in order to fully reproduce the shape of the data is also discussed. 
The DVCS cross sections in this paper represent the final set of experimental results from EOO-110, 
superseding the previous publication. 
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I. INTRODUCTION 

In the past two decades the Deeply Virtual 
Compton Scattering (DVCS) and Deep Vir¬ 
tual Meson Production (DVMP) reactions have 
emerged as powerful probes of the quark-gluon 
structure of the proton, neutron, and other 
atomic nuclei. DVCS, specifically, refers to the 
reaction "f*p —>■ py in the Bjorken limit of Deep 
Inelastic Scattering (DIS) j*p kinematics, but 
at low net invariant momentum transfer t to the 
target (in this case the proton). Experimentally, 
we can access DVCS through exclusive electro¬ 
production of real photons ep epy, where the 
DVCS amplitude interferes with the so-called 
Bethe-Heitler (BH) process (Figure]^. Given 
the knowledge of the proton elastic form factors, 
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the BH contribution is calculable in Quantum 
Electrodynamics (QED) since it corresponds to 
the emission of the photon by the incoming or 
the outgoing electron. 



FIG. 1. Lowest-order QED amplitude for the 
ep —^ ep 7 reaction including its decomposition. The 
momentum four-vectors of all external particles are 
labeled at left. The net four-momentum transfer 
to the proton is = (q — q')^ = (p' — p)^. In 
the virtual Compton scattering (VCS) amplitude, 
the (spacelike) virtuality of the incident photon is 
= —q^ = —(k — k')^ . In the Bethe-Heitler 
(BH) amplitude, the virtuality of the incident pho¬ 
ton is — = —t. Standard (e, e') invariants are 

Se = (k + p)'^, xb = Q^l{2q-p) and W'^ = (q + p)'^. 


DVCS is the simplest probe of a new class 
of light-cone matrix elements, called General¬ 
ized Parton Distributions (GPDs) [1]. These 
reactions offer the exciting prospect to obtain 
3-dimensional tomographic images of the trans¬ 
verse spatial distributions of partons (elemen¬ 
tary quarks and gluons) as functions of the par- 
ton light-cone momentum fraction m- In the 
kinematics of the present experiment, the GPDs 
are dominated by the quark light-cone matrix 
elements. The correlation of transverse spatial 
and longitudinal momentum information con¬ 
tained in the GPDs provides a new tool to eval¬ 
uate the contribution of quark orbital angular 
momentum to the proton spin [3] . 

The proof of Quantum Ghromodynamics 
(QGD) factorization theorems [ZlIH] established 
that in the Bjorken limit of high at fixed xb, 
the GPDs are the leading contribution to the 
7 *p —>■ 7 p amplitude in an expansion in inverse 
powers of (twist-expansion). Higher terms 
in the expansion are sensitive to more compli¬ 
cated correlation matrix elements {e.g. qqg cor¬ 
relations). The leading order DVGS amplitude 
is determined by four GPDs, which are defined 
in terms of vector H and E and axial vector 




FIG. 2. The leading-order Virtual Compton Scat¬ 
tering (VCS) amplitude in the limit of large , 
fixed Xb, and small t — A^. The kinematic variable 
^ = -{q + q'f /{2{q + q') ■ P), with P = {p + p')/2. 
In the aforementioned limit ^ —>■ xbI(2 — Xb) and 
2C ^ A+/P+ = (A° -b AQ/(P° -bPQ, with the z- 
direction parallel to P in the q -f P = 0 CM frame. 
Similarly, in the middle diagram, the quark and pro¬ 
ton lines are labeled by their ‘-I-’ momentum frac¬ 
tions and ’-b’ momentum components, respectively 


H and E light-cone matrix elements. The Gen¬ 
eralized Parton Distributions enter the DVGS 
cross section through Gompton Form Factors 
(GFFs), which are integrals over the quark loops 
of the two diagrams of Fig. For example, the 
GFF H corresponding to the GPD H is defined 
through (/ G {u,d, s}) 0: 


=E S i)] 



1 

^-x 


1 

C + X 


Hf{x,^,t) 


( 1 ) 


Thus, the imaginary part accesses GPDs along 
the line x = ±^, whereas the real part probes 
GPD integrals over x. The ‘diagonal’ GPD, 
is not a positive-definite prob¬ 
ability density, however it is a transition den¬ 
sity with the momentum transfer Aj_ Fourier- 
conjugate to the transverse distance r between 
the active parton and the center-of-momentum 
of the spectator partons in the target [TD]. Fur¬ 
thermore, the real part of the Gompton Form 
Factor is determined by a dispersion integral 
over the diagonal x = plus a D-teiui in- 
ill]. This D-teim |15j only has support in the 
ERBL region |a:| < ^ in which the GPD is de¬ 
termined by qq exchange in the t-channel. 
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GPDs have generated an intense experimen¬ 
tal activity. Beam spin asymmetries for DVCS 
in the valence region were first measured by 
the HERMES [TB] and CLAS [T7] collabora¬ 
tions. Cross sections were first measured at 
low-a:^ by the HI [18] and ZEUS collabora¬ 
tions m- These results were followed with 
more detailed studies of the Q^-, tU^-, and t- 
dependence of the cross sections The 

HERMES collaboration has measured a diverse 
range of asymmetries on the proton, including 
longitudinal-spin [21] , transverse-spin [5^ [2S] , 
beam-charge E3, and kinematically complete 
beam-spin asymmetries [28]. Detailed studies 
of the DVCS cross section as a function of 
<5^, and t by the ZEUS [TS] and HI Collabo¬ 
rations [21] demonstrated the factorization of 
the cross section, and the dominance of gluon 
CPDs at low xb- The first measurements of 
the DVCS cross section in the valence region 
were obtained by the present experiment |29j . 
together with an extraction of DVCS off the 
neutron [5^. A subsequent JLab Hall A exper¬ 
iment is analyzing the beam energy dependence 
of the DVCS cross section [3I1I32]. Beam-spin 
and longitudinal target spin asymmetries in the 
valence region were measured in CLAS [55H5S] . 
Extensive DVCS data taking has now started 
with the JLab upgrade m- Over the next few 
years, a broad CPD program is planned for the 
12 CeV beams at Jefferson Lab and the high en¬ 
ergy muon beams at CERN in the COMPASS 
experiment. 


II. THEORETICAL FRAMEWORK 


The photon electroproduction cross section of 
a polarized lepton beam of energy k off an unpo¬ 
larized target of mass M is sensitive to the co¬ 
herent interference of the DVCS amplitude with 
the Bethe-Heitler amplitude (see Fig. [^. It can 
be written as: 


dPa{X, ±e) 
dQ'^dx B dtd4>d(j)e 

\'j-bh 


d'^ao 1 
dQ'^dxB e® 

|2 \^DVCSl'^ 




(2) 


_ OfQED 1 

dQ'^dxB 167r2(se — M‘^)‘^xb Vl + 

=4M2a;|/Q2 

Se = 2Mk + 


where 4>f. is the azimuthal angle of the scattered 
electron around the beam axis in the labora¬ 
tory frame, (j) is the azimuthal angle between 
the leptonic and hadronic planes defined in the 
Trento convention [35] , A is the electron helicity 
and the -l-(—) stands for the sign of the charge 
of the lepton beam. The cross section does not 
depend on ipf, and this angle is integrated over, 
leaving effectively a 4-difFerential cross section. 
The BH contribution is calculable in QED, 
given the ^1% knowledge of the proton elastic 
form factors in our range of — < < 0.4 CeV^. The 
other two contributions to the cross section, the 
interference term I and the DVCS squared term 
l-y-DVCSl ^ provide complementary information 
on CPDs. It is possible to exploit the structure 
of the cross section as a function of the angle 
(p to separate up to a certain degree the differ¬ 
ent contributions to the total cross section [33] . 
The BH term is given in [9], Eq. (25), and only 
its general form is reproduced here: 


SLo Cn ^ cos(n((') 

xlty^l + e^yVi{cP)V2{cP) 


The harmonic terms depend upon bilin¬ 
ear combinations of the ordinary elastic form 
factors Fi{t) and F 2 {t) of the proton. The fac¬ 
tors Vi are the electron propagators in the BH 
amplitude [3]. 

The interference term in Eq. ([^ is a linear 
combination of GPDs, whereas the DVCS^ term 
is a bilinear combination of GPDs. These terms 
have the following harmonic structure: 
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1 = 


XBy^'Pi{4>)V2{(j))t 


4 + cos(n(/)) + Xsl sin(n</))] 


(5) 


n—1 


\T 


■Dyes 12 


nves 


y'^Q'^ 


2 

E 

n—1 


^DVeS 


cos(n^) + sin(n0)] 


( 6 ) 


The and (c, s)^ harmonics are domi¬ 

nated by twist-two GPD terms, although they 
do have higher-twist admixtures that must be 
quantified by the dependence of each har¬ 
monic. The (c, and (c, s)f harmon¬ 

ics are dominated by twist-three matrix ele¬ 
ments, although the same twist-two GPD terms 
also contribute (but with smaller kinematic 
coefficients than in the lower Fourier terms). 
The (c, 3 ) 2 ^'^^ and (c, s)f harmonics stem 
only from twist-two double helicity-flip gluonic 
GPDs. They are formally suppressed by Og and 
will be neglected here, but they do not mix with 
the twist-two quark amplitudes. 

The bilinear DVGS term has a twist-2 contri¬ 
bution that reads: 


„Dycs _ 2 ^ ‘^y + y^ + 2 y^ 


1-f 


^UVCS(jr,jF*) 

where T represents the set {H, £,'H, £} of 
twist-2 GFFs. The Fourier coefficients and 
of the interference term are: 

cj = C^+3fieC^’”(J-) + C'o\3fie44”(J-eff), 


+ 5o\3m5o^4"(J-eff). 




.X,n / 


( 8 ) 


The above coefficients are defined in terms of 
the photon helicity-conserving 

^V,n ^A,n 


-I-+ 


c 


-I-+ 


(9) 


QV,n qA^u 

{T)=C\F) + 


( 10 ) 


r 


and helicity-changing amplitudes 


CoT(-^eff) = 


V 2 K 


2 — Xb Q 


C^(J'eff) + 


a 


V,n 
0+ r^xy 


a 




(Teff) + 


0 + 






0 + 


( 11 ) 


V2 K 

2 — Xb Q 


C^(^eff) + 

A^n 


-yV,n 


•"0+ 


^1 


0 + 


( 12 ) 


The complete expressions of kinematic coeffi¬ 
cients C"j, S^i, and K are given in [iD]. The 
and terms are respectively linear and 

bilinear combination of GFFs. For example: 

c^(^:F) = F,n+aFi+F2)n-^F2£. ( 13 ) 


III. THE EOO-110 EXPERIMENT 

The EOO-110 [H] experiment ran in Hall A at 
Jefferson Lab in the fall of 2004. Its goal was to 
measure the dependence of the DVGS [42] 
helicity-dependent cross sections at fixed value 
of Xb: 


(ra = 


A^cr = 


d‘^a{X = -1-1) ^ d‘^(7{X = —1) 
dQ'^dxBdtdcf) dQ'^dxBdtdcf) 
d^a{X = -fl) d'^a{X = —1) 


dQ'^dxBdtdcj) dQ'^dxBdtdcf) 


(14) 


(15) 


Tab. U summarizes all 
tings of this experiment. 


the kinematic set- 
in addition to the 
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FIG. 3. (Color online) Distribution of H{e,e'y)X 
events in the [xb, Q^] plane, for Kin2 (xb = 
0.36, = 1.9 GeV^) and Kin3 {xb = 0.36, 

= 2.3 GeV^). Events for KinX2 {xb = 0.39, 
= 2.06 GeV®) and KinX3 {xb = 0.34, = 

2.17 GeV^) are bounded by the two horizontal lines 
at = 1.95 GeV^ and = 2.30 GeV^ 


dependence at fixed xb, we present here 
new results on the x^-dependence of the DVCS 
cross section at fixed by using a subset 
of the data from the Kin2 and Kin3 settings 
with 1.95 < < 2.30 GeV^, as illustrated in 

Fig. [31 We labelled these new settings Kin X2 
and Kin X3. They are centered at xb = 0.40 
and Xb = 0.34, respectively, for an averaged 
= 2.1 GeV2. 

The setup of the experiment is shown in 
Fig. i In order to counter the small cross 
section, this experiment used the high lumi¬ 
nosity in Hall A of Jefferson Lab, running at 
10^^ cm“^s“^, which corresponds to 2.25 /rA of 
electron beam on a 15-cm-long liquid hydrogen 
target. The scattered electron was detected in 
the Hall A left High Resolution Spectrometer 
(HRS), which provides a momentum resolution 
5p/p = 2 • 10“^ and an angular resolution of 
2 mrad in the horizontal plane |43j . This pin¬ 
points the electron kinematics {xb and Q^), the 
electron scattering plane, and the momentum 
direction q of the virtual photon of the virtual 
Compton amplitude. The emitted photon was 
detected in an electromagnetic calorimeter cov¬ 
ering ^ 0.1 sr, with its front face 1.1m from the 
target center, and centered in the direction of 


the virtual photon (shifted by half a calorime¬ 
ter block). The spectrometer acceptance of 6 
msr and ±4.5% in momentum selects virtual 
photons in a small solid angle or ~ 3 msr, as 
illustrated in Fig. The detected photon di¬ 
rection (two angles) with respect to the virtual 
photon direction (calculated using the electron 
kinematics) determines the remaining two kine¬ 
matic variables of the reaction: t and </>. The 
measurement of the detected photon energy al¬ 
lows for an exclusivity cut based on the squared 
missing mass of the recoil proton. As a cross¬ 
check on exclusivity, the recoil proton was de¬ 
tected in the Proton Array, a set of 100 blocks 
of plastic scintillator in a C-ring configuration 
around the virtual photon direction. This ge¬ 
ometry was selected in order to have a simple 
azimuthal symmetry around the virtual photon 
direction, which is a key element for a smooth 
(j) acceptance. 

The basic equipment of Hall A, including the 
beamline, target system, and dual spectrome¬ 
ters is described in [43]. The following sections 
provide details specific to the present experi¬ 
ment. 


HRS 



FIG. 4. Setup of the EOO-110 experiment in Hall A 
of Jefferson Lab. The photon calorimeter as well 
as the proton array were centered on the virtual 
photon direction q, then shifted sidewise by half a 
calorimeter block away from the beam to limit the 
singles rate on the detector elements close to the 
beamline. 
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Setting 

k' (GeV/c) 

e. n 

(GeV"*) 

Xb 

0, n 

W (GeV) 

E-, (GeV) 

Kinl 

3.53 

15.6 

1.5 

0.36 

-22.3 

1.9 

2.14 

Kin2 

2.94 

19.3 

1.9 

0.36 

-18.3 

2.0 

2.73 

Kin3 

2.34 

23.8 

2.3 

0.36 

-14.8 

2.2 

3.32 

KinX2 

2.94 

20.1 

2.06 

0.39 

-18.6 

2.03 

2.71 

KinX3 

2.36 

23.1 

2.17 

0.34 

-14.5 

2.26 

3.33 


TABLE I. Experimental ep — >■ ep 7 kinematics, for incident beam energy Et, = 5.7572 GeV. 9q is the central 
value of the q- vector direction. E-^ is the photon energy for t = tmin. A subset of Kin2 and Kin3, with the 
cuts shown in Fig. provides the kinematic settings KinX2 and KinX3 at fixed but varying xb- Note 
that only the average kinematics for each setting are listed in this table : in order to minimize systematic 
bin centering effects, the results are presented or listed using the kinematics of each bin in xb, and t 
according to their averaged experimental value in the bin. Our extraction procedure ensures that all ij> bins 
are evaluated at the same kinematic setting, as explained in section 


A. Electron Beam 

1. Beam Energy 

The incident beam energy is measured by de¬ 
termining its bend in the arc section of the Hall 
A beamline [H]. Its deflection angle is com¬ 
puted from a set of wire scanners. The mag¬ 
netic field integral of the eight dipoles of the 
beamline is compared to a reference magnet 
(9*^ dipole). The measurement of the beam en¬ 
ergy made during the experiment resulted in the 
value Et = 5757.2 ± O.lgtat ± O.lgyst MeV. 


2. Beam Current 

The beam current is measured using two res¬ 
onant RF cavity monitors (Beam Current Mon¬ 
itors) tuned at the frequency of the accelerator 
(1.497GHz). The voltage at their outputs is 
proportional to the beam current and provides 
a continous monitoring of its value during the 
experiment. The absolute reference is provided 
by a separate monitor, a Parametric Current 
Transformer [43] , which is calibrated by passing 
current of known value through a wire inside the 
beam pipe. 


3. Beam Polarization 

The electron beam polarization was measured 
concurrently with the regular data taking using 
the Hall A Compton polarimeter [45]. At the 
entrance of the Hall, the beam is deflected by a 
chicane and interacts with a circularly polarized 
photon beam. The polarization of the electron 
beam can be obtained from the counting rate 
asymmetry from opposite beam helicities. The 
electrons that interact with the photon beam 
are detected by silicon micro-strips, while those 
that do not interact continue towards the exper¬ 
imental target. The photon beam is provided 
by a resonant Fabry-Perot cavity that ampli¬ 
fies a 230 mW Nd:YaG laser (A = 1064 nm) to 
1200 W. 

The statistical error of a Compton measure¬ 
ment is inversely proportional to the square root 
of the number of events and to the analyzing 
power of the polarimeter. In this experiment, 
a 1% statistical error could be achieved in 2.5 h 
of data taking. However, this is far from being 
the limiting factor. Since the Compton data was 
taken during normal DVCS running, we can av¬ 
erage over long periods of time in order to make 
the statistical error negligible. 

Beam polarization results can be readily ob¬ 
tained from the electron detector. The electron 
detector consists of 4 planes of 48 silicon micro¬ 
strips, standing 4.6 mm above the beam axis 













FIG. 5. (Color online) Perspective view of the 
downstream face of electromagnetic calorimeter. 
The virtual photon ( 7 *) acceptance is shown pro¬ 
jected on the calorimeter plane. The distance (black 
line) between the ’’impact” of the virtual photon 
and the detected real photon ( 7 ) position is roughly 
proportional to \/t — tmin for the small tmin — t val¬ 
ues of this experiment. The angle <j) is the azimuthal 
angle of the photon with respect to the plane formed 
by the incident beam and the virtual photon. For 
the central (e, e') kinematics indicated by 7 *, (j) is 
measured counter-clockwise from the horizontal (to 
the right) direction. 


during the DVCS experiment. Figure § (top) 
shows the electron counting rate versus strip 
number in one of the detector planes for a typi¬ 
cal Compton run of 3 h duration. The detector 
is located behind the third dipole of the Comp¬ 
ton chicane and the strip number gives the po¬ 
sition of the scattered electron along the dis¬ 
persive axis with a resolution of 200 /rm. Hence 
the horizontal axis of the plot is proportional 
to the energy lost by the electron (and given 
to the photon). The Compton energy spec¬ 
trum shows up as a flat rate on the first strips. 
The background spectrum has a 1 /E shape, like 
bremsstrahlung. The differential asymmetry as 
a function of the electron energy (strip number) 
is shown in the bottom plot of Fig. for each 
of the laser polarization states. 

The systematic error in the polarization mea¬ 


surement due to the uncertainty of the laser 
polarization is 0.7%. The maximum devia¬ 
tion of scattered electrons for a beam energy 
of ~5.75GeV is 21.5 mm at the electron detec¬ 
tor plane, which makes a calibration error of 
200/rm/21.5mm=0.93%, which propagates to 
1.9% to the polarization measurement. The to¬ 
tal systematic error associated to the beam po¬ 
larization measurement is 2%. 

Figure shows the Compton polarimeter re¬ 
sults for the full experiment duration, where 
only the electron detector was used in the anal¬ 
ysis. The beam polarization was 75.3±0.1stat i 
2.0syst% in average during the experiment. The 
low polarization values at the beginning of the 
experiment (first 3 points in Fig. correspond 
to the period when the polarization was not yet 
optimized for Hall A. 


B. Liquid Hydrogen Target 

The standard Hall A cryogenic target sys¬ 
tem |46j was mounted inside a scattering cham¬ 
ber custom built for the EOO-110 experiment. 
Its walls were thinner than those of the usual 
Hall A scattering chamber. Also, a larger exit 
beam pipe was constructed. The scattering 
chamber was made of a 1 cm spherical shell of 
aluminum, allowing for low energy protons to go 
through (minimum momentum of 305 MeV/c, 
corresponding to a cut on the kinematic variable 
t of -0.091 GeV^). Moreover, the new scattering 
chamber accommodates the spherical symmetry 
of the reaction and makes energy losses nearly 
independent of the scattering angle, except for 
extended target effects. The larger exit beam 
pipe reduces the background. 

The cryogenic target has three target loops, 
two of which were used for the DVCS experi¬ 
ment: a liquid hydrogen (LH 2 ) loop and a liq¬ 
uid deuterium (LD 2 ) loop [17]. Each of the two 
liquid loops had an aluminum cylindrical target 
cell, 15 cm long. An additional solid target lad¬ 
der was attached to the system for calibration 
purposes. The targets are arranged in a verti¬ 
cal stack, which can be moved from one position 
to another by remote control. The solid target 
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FIG. 6. (Color online) Signal and background rates 
(normalized to the electron beam current) in one of 
the planes of the electron detector as a function of 
strip number (top), and asymmetry measured for 
each of the laser polarization states and the back¬ 
ground as a function of strip number (bottom). 
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{P,} = 75.32% ± 0.07% 


{S/N) = 27.60 ±4.12 
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FIG. 7. (Color online) Compton polarimeter re¬ 
sults, using only the electron detector. Beam po¬ 
larization is shown in the upper plot. Lower plot 
shows the signal to background ratio. The first 
three points in the beginning of the experiment cor¬ 
respond to a non-optimal Wien angle setting. 


• C: 1 mm thick carbon, serving as a point¬ 
like target. 

• Empty: Position used to reduce radia¬ 
tion on detectors while beam was used for 
other purposes (beam size measurements 
using wire scanners and other beam tun¬ 
ings). 


C. Hall A Spectrometer 


ladder contained the following targets: 

• Optics: Seven 1-mm-thick carbon foils 
used for optics calibration of the HRS. 

• Two dummy targets: ±2 cm and ±7.5 cm 
A1 foils to study target walls effects. 

• Cross hair: Aluminum foil with a milled 
cross, used to measure beam position with 
respect to the target. 

• BeO: It allows to see the beam spot at 
the target through a camera installed in 
the scattering chamber. 


The High Resolution Spectrometers (HRS) in 
Hall A consist of 4 superconducting magnets in 
the configuration QQDQ. In the EOO-IIO ex¬ 
periment, the left HRS was used to detect the 
scattered electron and therefore define the vir¬ 
tual photon kinematics in an accurate way. The 
main components of the detector stack are as 
follows: a set of two scintillator planes called 
SI and S2m giving very fast and good timing 
signals; two vertical drift chambers for track re¬ 
construction; a gas Cerenkov counter for n/e 
discrimination and a Pion Rejector composed 
of two layers of Lead Glass blocks which is used 
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in addition to the Cerenkov detector to select 
a clean sample of electrons. A fast signal from 
S2m in coincidence with the Cerenkov detec¬ 
tor was used as a level 1 trigger for the rest of 
the electronics. It is useful to recall the angu¬ 
lar acceptance of the left HRS for electrons : 
±30 mrad horizontal, ±60 mrad vertical, and 
±4.5% in momentum. 


D. Calorimeter 

One of the key elements of this experiment 
was a dedicated electromagnetic calorimeter, 
consisting of a 11 x 12 array of lead fluoride crys¬ 
tals, each 3x3xl8.6cm^. The crystals were 
purchased from SICCAS (Shanghai). PbF 2 
was selected as a pure Cerenkov medium, to 
minimize hadronic backgrounds and to obtain 
the shortest possible signal without exponen¬ 
tial tails. The size of the blocks is adapted 
to the radiation length and Moliere radius of 
PbF 2 so that a shower is almost completely 
contained in a cluster of 9 blocks, both lon¬ 
gitudinally and transversally. Each block was 
equipped with a Hamamatsu R7700 fine-mesh 
photomultiplier tube (PMT). During the exper¬ 
iment, the relative gains of the PMTs were peri¬ 
odically monitored using a cluster of LEDs that 
could be moved across the calorimeter face on 
an X-Y stage. However, the large luminosity of 
the experiment induced radiation damage near 
the front face of the blocks. Since the LED 
light was injected in the front face, the LED 
method proved unreliable to measure the true 
signal variation of high energy photon or elec¬ 
tron showers. Indeed, the Cerenkov light from 
a multi-GeV y-ray is mostly produced deeper 
in the crystal, avoiding most of the damaged 
area. Figure]^ shows simulated shower profiles 
for 4 GeV electrons, at a typical energy for the 
electromagnetic background and for photons of 
various energies expected from DVCS photons. 

The crystal-by-crystal calibration coeffi¬ 
cients were obtained from kinematically over¬ 
constrained elastic scattering: H(e, 
in which the electron is detected in the calorime¬ 



FIG. 8. (Color online) Shower longitudinal shower 
profile for different incoming particles into the 
DVCS calorimeter, obtained from the Monte-Carlo 
simulation. The aluminum shielding corresponds to 
the total thickness of material between the target 
and the calorimeter crystals, and includes both the 
scattering chamber and some additional A1 shield¬ 
ing in front of the calorimeter front face. 


ter and the proton is detected in the HRS. In 
order to illuminate the full acceptance of the 
calorimeter with elastic electrons, it was neces¬ 
sary to move the calorimeter back to a distance 
of 5.5 m from the target center during these 
runs. Data at 1 m, covering only the center 
part of the calorimeter were taken additionally 
as a consistency check. Figure shows the en¬ 
ergy resolution of the calorimeter as measured 
during the elastic calibration runs. Two elastic 
calibrations were made, one a few weeks after 
the start of the experiment and another one a 
few weeks before it finished. The cafibration 
coefficients changed by a considerable amount 
for some blocks, but the energy resolution did 
not degrade during the almost 3 months of data 
taking. Since calibration coefficients changed 
with time, in order to keep a good energy reso¬ 
lution all along the experiment, we interpolated 
these coefficients between the two cafibrations 
runs, together with an extrapolation before and 
after them. This was done based on the ra¬ 
diation dose accumulated by each block. This 
dose is proportional to the beam current and de¬ 
pends on the block polar angle with respect to 
the beam line and also on the target type (LH 2 
or LD 2 ). The relative dose accumulation for 
each block was estimated from its PMT anode 
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FIG. 9. (Color online) Energy measured in the 
calorimeter minus energy expected from elastic 
kinematics during elastic calibrations runs. In both 
elastic calibration periods, we obtained 2.4% energy 
resolution at an elastic energy of 4.2 GeV. The re¬ 
sults of the second calibration when first calibration 
coefficients are used are also plotted to show the ne¬ 
cessity of a careful monitoring of the coefficients in 
between these two calibration points. 


current monitoring |48j . In addition, the cali¬ 
brations were monitored in situ using the miss¬ 
ing mass peak of the reaction D(e, eQaio^HRs)TT 
and both the missing mass and ( 77 ) invariant 
mass peaks of the reaction H(e, e' 7 r°)p. Overall, 
the calibration coefficients were known for any 
given time at the 1 %-level. 


E. Proton Array 

In order to detect the full exclusive final state, 
a recoil detector was built to tag the DVCS pro¬ 
ton. The recoil proton direction for an exclusive 
event can be inferred from the information of 
the HRS and the calorimeter, therefore one can 
check in the Proton Array (PA) if the proton 
was actually at the right position. The main 
difficulties of such a detector is that it needs 
to detect low-momentum protons in a large ac¬ 
ceptance, close to the beam line, with as high 
an efficiency as possible. The PA subtended an 
acceptance (relative to the nominal direction of 
the virtual photon) of 18° < < 38° and 


45° < 4’-y*p < 315°, arranged in 5 rings of 20 
detectors as shown in Fig. The scintillator 
blocks were fabricated by Eljen Technology as 
5 distinct tapered trapezoids, each 30 cm-long, 
in order to form a hermetic ring pointing at the 
target center. Each scintillator is equipped with 
a Photonis XP2972 PMT and a custom voltage 
divider/pre-amplifier circuit. This allowed us 
to operate the PMTs at low gain, to accommo¬ 
date the high backgrounds in this open geome¬ 
try. The 90° cut-off in (j)j*p corresponds to the 
exit-beam pipe in the kinematic setting where 
the detector stack is the closest to the beamline. 


F. Sampling Electronics 

The EOO-110 experiment was designed with 
open detectors at low angles (the blocks of the 
calorimeter closest to the beam line were at 
6.5°) with limited shielding running at high lu¬ 
minosity. High singles rates up to 10 MHz 
were expected and also measured in a test 
run during the design phase. In this environ¬ 
ment, regular ADCs even with a reduced gate 
are strongly affected by pile-up. We therefore 
chose to use digitizing electronics for all the 
electronic channels of the dedicated detectors 
(PbF 2 calorimeter and Proton Array), namely 
a custom 6 U 16 VME (A24/D32) module sam¬ 
pling system based on the Analog Ring Sampler 
(ARS) CMOS ASIC developed at CEA-Saclay 

mm- 

The ARS uses the concept of analog mem¬ 
ories to sample data at a clock rate of 1 GHz: 
each channel contains a circular array of 128 ca¬ 
pacitors: every 1 ns, the ARS points the signal 
to the next capacitor, eventually overwriting it¬ 
self after 128 ns. When a trigger is issued, the 
capacitor array is isolated and the previous 128 
samples are stored. During the next 500 ns, a 
separate trigger module (described below) de¬ 
cides whether or not to digitize the event, hol¬ 
lowing a validation from the trigger, each capac¬ 
itor array is digitized in parallel using a 12 -bit 
Flash ADC at a rate of 1 /rs per sample, for a to¬ 
tal of 128 fis per channel. During this long dig¬ 
itization period, we observe an exponential RC 
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decay of the samples. This is compensated by a 
stable baseline included in the pulse waveform 
analysis. Each ARS ASIC contains four chan¬ 
nels, and four ARS ASICs were implemented 
onto each VME board for a total of 16 channels 
per board. Figure 10 shows a typical calorime¬ 
ter signal as a function of time, read out by the 
ARS system. 



FIG. 10. Flash ADC value as a function of time 
recorded by the ARS system for a typical calorime¬ 
ter pulse. 


G. Trigger 

The Data Acquisition (DAQ) trigger for this 
experiment was a two-level system. A standard 
HRS electron trigger was formed from the coin¬ 
cidence of the Cerenkov and S2m signals. This 
Level-1 signal generated the “Stop” to freeze the 
analog data in the ARS. The “Validation”, or 
Level-2 signal is generated by a dedicated DVCS 
Trigger module. The DVCS Trigger includes a 
large backplane, containing FPGA logic. Each 
PbF 2 signal is first sent to a Trigger Daughter 
card, where it is split with one branch going to 
an ARS input and the second branch passing to 
a Fast ADC chip on the Daughter card. Each 
Daughter card has 4 channels. The 132 trig¬ 
ger ADCs are gated by the Level-1 signal, with 
a programmable width generally set to 60 ns. 
Following digitization, the FPGA logic forms lo¬ 
cal 2x2 overlapping cluster sums. If a cluster 
sum is found above a programmable threshold 


(generally set to 1 GeV equivalent), then the 
Level-2 Validation signal is set to true. This is 
completed within ~ 500 ns. In the absence of 
a validation signal at the end the 500 ns win¬ 
dow, a fast clear is issued to the DVCS Trig¬ 
ger and ARS. In this way, for random Level-1 
triggers, the deadtime is only ^ 500 ns, and 
the full readout is incurred only for a genuine 
H(e, e' 7 )X coincidence (including accidentals). 
Notice that only calorimeter channels belonging 
to 2 X 2 clusters above threshold were digitized 
and recorded for each event. 


The PA was not in the trigger and was read 
at every HRS-calorimeter coincidence. As the 
virtual photon has an almost fixed direction, 
the approximate region of the proton detector 
that a DVCS proton would hit, can be inferred 
from only the calorimeter region hit by the pho¬ 
ton. The list of proton detector channels to 
be read out can then be made on the fly out 
of the information provided by the calorimeter 
trigger module. The accuracy of this predic¬ 
tion is dominated by a convolution of the HRS 
acceptance and the calorimeter energy resolu¬ 
tion. Processes such as multiple scattering in 
the target and the scattering chamber can also 
affect the accuracy of the prediction. The look¬ 
up table of proton detector channels to read as 
a function of calorimeter ones was computed us¬ 
ing the Monte-Carlo simulation, which included 
a realistic description of all the elements of the 
experimental setup (see section IVE). The de¬ 
tector inefficiency for DVCS protons due to this 
online readout choice is smaller than 0.1%. De¬ 
pending on the kinematic setting, only between 
15% and 30% of the proton detector blocks 
needed to be read out in average. This re¬ 
duced the amount of data to record and there¬ 
fore the acquisition deadtime. Note that no 
threshold was set in any of the proton detector 
channels, so that even very low-energy protons 
could be detected. The communication between 
the calorimeter and proton array crates neces¬ 
sary for this block selection in the proton detec¬ 
tor was made possible by a custom multiplexer 
module (MUX) which allowed the calorimeter 
trigger module to send its data to the proton 
array crate. 
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The HRS DAQ functions in the standard way 
for a Hall A coincidence experiment, with all 
HRS analog PMT signals sent through delay 
cables corresponding to 880 ns. Scintillator sig¬ 
nals go in common start into a LeCroy 1875 
high resolution TDC. The delayed HRS signals 
arrive at their respective ADC and TDC inputs 
after the Level-2 decision is made. 

Even though the ARS represent a consider¬ 
able advantage for this type of experiment it 
has an obvious drawback: the amount of data 
to transfer is about a hundred times higher 
than a regular HRS event. If the 232 ARS 
channels of the EOO-110 experiment were to be 
recorded at every event, the event size would 
be 232x 128x 16/8~60kB. Typically, only ~ 40 
ARS channels were recorded at every event. 


IV. DATA ANALYSIS 

As mentioned before, the selection of the 
ep —)■ ep 7 final state is based on a missing- 
mass analysis of the ep —> ejX event sample. 
This is made possible by the excellent momen¬ 
tum resolution of the Hall A HRS and the fair 
energy and position resolutions of our dedicated 
electromagnetic calorimeter. The following sub¬ 
sections describe the selection of electron can¬ 
didates from the HRS, and the analysis of the 
calorimeter in order to select the final-state pho¬ 
ton. We will then focus on the final steps to en¬ 
sure that our ep ep'y selection is efficient and 
its purity close to perfect. Finally, we will de¬ 
scribe the normalization procedure, the Monte- 
Carlo simulation and the method used to ex¬ 
tract cross sections from our data and Monte- 
Carlo events. 


A. HRS Analysis 


The HRS Cerenkov detector was used for the 
electron identification. The number of photo¬ 
electrons detected is 7 on average so that the 
distribution is Poissonian. Figure 11 shows the 
distribution of the sum of all 10 PMTs (in ADC 
channels). The first ’peak’ in the spectrum is 



FIG. 11. (Color online) Distribution of the sum of 
all 10 Cerenkov PMT ADC values, for each kine¬ 
matic setting. The cut applied to remove the 1- 
photoelectron signal from data is also shown. 


the tail of the electronic noise in the pedestal. 
We remove 1-photoelectron events (either ther¬ 
mal emission in the PMT or 5-rays from pi- 
ons) by applying a cut at 150 ADC channels. 
The 1-photoelectron peak is only visible if the 
Cerenkov signal is removed from the trigger, 
and a cut is made on the Pion Rejector to select 
minimum ionizing particles {i.e. pions). 

Figure shows the distribution of the re¬ 
action point along the beam Vz reconstructed 
by the HRS. The target center relative to the 
Hall center was determined to be 7.8 mm down¬ 
stream. A cut in order to avoid the contribution 
from the target cell wall was applied to the data: 
—6.00 cm < Vz < 7.50 cm. 

Figure shows the resolution on the ver¬ 
tex reconstruction as measured with a carbon 
multi-foil target. The thickness of each foil is 
1 mm, and the HRS was at 37.69° during this 
run. The HRS vertex resolution varies as: 


sin6*HRS ' 


(16) 


The a measured at 37.69° is 1.87 mm, which 
means 1.2 mm at 90° (the a introduced by the 
foil thickness is (l/vT2)mm and can be ne¬ 
glected) . 

The HRS acceptance is a hypervolume de¬ 
pending on 5 correlated variables: Xtg and 9tg 
(the position of the particle and the tangent of 
the angle made by its trajectory along the dis- 







14 



FIG. 12. (Color online) Reaction point along the 
beam reconstrncted by the HRS. The cut on the 
target length applied is shown by the vertical lines. 
The 7.8 mm downstream shift of the target observed 
during the experiment is also evident. 



FIG. 13. (Color online) Top: Resolution of the ver¬ 
tex reconstruction from a multi-foil target. Bottom: 
Close-up of the central foil fit, which results in a 
a = 1.9 mm resolntion. The foil thickness is 1mm 
and the HRS was at 37.69° dnring this run. 


section and higher-momentum particles tending 
to fly closer to its shorter base (high magnetic 
field) side, which makes the accepted range of 
smaller for higher btg. Finally, increasing ytg 
requires decreasing (j)tg in order for the particle 
to get into the spectrometer entrance window. 
Making cuts independently in each of the vari¬ 
ables to limit events to flat acceptance regions 
in each of them is thus very inefficient. Instead, 
we used an acceptance function m, which al¬ 
lows to place a 4-dimensional cut {xtg = 0 is 
assumed). This procedure is almost twice more 
efficient than the traditional sequential accep¬ 
tance cuts. This function takes the arguments 
Vtg, (^tgi 'Ptg and 5tg and returns a so-called R- 
value which is the minimum distance (in radi¬ 
ans) to the {9tg,4’tg) solid angle acceptance re¬ 
gion appropriate for a given value of ytg and 
5tg. A value of 5 mrad was used in order to 
constrain a well-defined region of the HRS ac¬ 
ceptance. The cross-section results varied by no 
more than 1% when increasing the R-function 
cut. This value was used as an estimate of the 
HRS acceptance systematic uncertainty. 


B. Calorimeter Analysis 

The calorimeter analysis is done in two steps: 
first, the recorded ARS waveforms are analyzed 
in order to extract the time and energy informa¬ 
tion. Then, an algorithm is used to aggregate 
the block information into photon clusters with 
a measured position, time and total energy. 


1. ARS Waveform Analysis 


persive direction), ytg and (j)tg (the position and 
the tangent of the angle perpendicular to the 
dispersive direction), and Stg (the fractional de¬ 
viation of the particle momentum with respect 
to the central momentum of the HRS). Trajec¬ 
tories of higher-momentum particles have lower 
curvature in the dipole, and in order for them to 
fit into the spectrometer they need to have lower 
6tg. The dipole magnet has a trapezoidal cross 


All the detector channels of the electromag¬ 
netic calorimeter were equipped with ARS elec¬ 
tronics, which allowed to save the full waveform 
of blocks that were recorded during a trigger, in 
a manner similar to a digital oscilloscope. In or¬ 
der to extract time and amplitude information 
from the ARS, a waveform analysis is needed 
which is performed offline. 

Each pulse as a function of time is described 
by a reference pulse multiplied by an amplitude. 
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For an ideal event without noise, the amplitude 
of the pulse and its arrival time are free parame¬ 
ters. For any given arrival time t, the amplitude 
a{t) which best fits the signal {xi} is simply 
given by the one which minimizes: 

imax 

X^it) = ^ (a;* - a{t)hi-t - b{t)f , (17) 

imin 

where {hi} is the reference shape. Notice that 
we also fit a flat baseline b{t). Reference shapes 
for each individual PMT are determined exper¬ 
imentally from data, using elastic calibration 
runs, where the probability of pile-up is very 
small. In order to reduce the impact of acci¬ 
dental events, only imax — *min=80 ARS sam¬ 
ples were used in the calorimeter analysis, cen¬ 
tered around the expected arrival time of DVCS 
events, which because of cable lengths, varies 
slightly from one channel to another. The par¬ 
tial derivatives of x^(t) with respect to a{t) and 
b(t) yield a linear set of equations in order to 
obtain the best amplitude for any given arrival 
time t. If the minimum value of x^(t) found 
for all the possible t is above a given analy¬ 
sis threshold yf, the algorithm will fit a second 
pulse to the waveform by minimizing: 

'i'max 

imin 

a2{ti,t2)h_t2 - b{ti,t2))‘^ , (18) 


The threshold value Xi used for the analysis 
corresponded to an effective missed pulse of 
~280 MeV for each particular calorimeter block 
(which translates to slightly different ARS chan¬ 
nel thresholds due to the different calibration 
of each block). Also, if the of a fit by a 
flat-line b was below an equivalent energy of 
Xo ~ 40 MeV, no pulse was fitted and the signal 
was discarded. Finally, if 2 pulses were found 
with a relative arrival time smaller than 4 ns, 
the algorithm returned the best single pulse fit 
since 2-pulse results proved to be unstable in 
those cases. 

The waveform analysis of the proton array 
ARS data used the same algorithm, but with 
slightly different parameters. Energy thresholds 
were set to yg ^ 2 MeV and y^ ~ 15 MeV in or¬ 
der to best ht the much smaller recoil proton 
energies in the detector. Due to the high count¬ 
ing rate in the detector, only 30 ARS samples 
were used for the fit. Also, time windows to 
search for pulses were set to —20 < ti, ^2 < 20 ns 
around the expected event signal. 

Overall, the waveform analysis of ARS signals 
increases the energy resolution in the DVCS 
calorimeter by a factor of 2-3 (depending on 
the background level) with respect to results ob¬ 
tained integrating the signal in a 60 ns window. 
We found about 8% of events in the calorimeter 
with some pile-up from accidentals. 

2. Clustering Algorithm 


for every combination of and ^2- For every 
pair of ti and O and the corresponding fitted 
amplitudes and baseline, a reduced y^ is also 
computed in a time window of ±20 ns around 
the minimum of the pulse. The minimum re¬ 
duced y^ found determines the amplitudes and 
arrival times of the pulses. Pulses were searched 
in a [—20, 25] ns interval around the expectation 
arrival time of events, in steps of I ns. An im¬ 
proved time resolution is obtained by interpo¬ 
lating around the time that minimizes the y^ 
for any time t = ti,t 2 '. 


The algorithm used to separate clusters in the 
electromagnetic calorimeter is based on a cellu¬ 
lar automata, as described in [52] , and uses only 
pulses arriving within a [-3,3] ns interval. This 
coincidence time window is more than 6 times 
the time resolution of the detector (~0.8ns). 
For each cluster found, the total photon energy 
E is taken to be the sum over the deposited 
energy Ei in each of the cluster blocks: 

E = Y,E^ E,= aA ,, ( 20 ) 


t = ti.Xmrn) + 


xi-1 


■Xt+1 


2 (y?+i ± Xt -1 - 2 yLn)' 


(19) 


where Ai is the signal amplitude collected in 
block i and Ci its calibration coefficient. The 
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impact position Xdus is calculated as the sum 
of blocks positions Xi weighted logarithmically 
by the relative energy deposition in each of them 

m- 


^clus — 


Wj 


Wi = max{0, Wq + h\[Ei/E)} . 

( 21 ) 


The parameter Wq allows a further tuning of 
the relative weight between blocks: as Wq oo 
the weighting becomes uniform regardless of the 
energy deposited in each block, whereas small 
values of Wq give a larger relative weight to 
blocks with large energy deposition. The value 
of Wq fixes the energy threshold for blocks to 
be taken into account in the position determi¬ 
nation: blocks with a relative energy deposition 
less than e~^° are neglected in the calculation. 

The calorimeter was placed at 110 cm from 
the 15-cm-long target. The incidence angle 
of particles on the front face of the calorime¬ 
ter could therefore vary by significant amounts: 
corrections due to the vertex position in the tar¬ 
get needed to be applied. Furthermore, the elec¬ 
tromagnetic shower does not begin at the sur¬ 
face of the calorimeter, but at a certain depth 
as shown in Fig. This depth is, to first ap¬ 
proximation, independent of the incident par¬ 
ticle energy. Taking these two effects into ac¬ 
count, the position Xdus given by equation (21) 
is corrected by: 


f ^ ® 

^corr — ^clus 1 1 / r 2 2 

V V ^vc + ^ 

where Lyc is the distance from the vertex to the 
calorimeter and a is the distance of the elec¬ 
tromagnetic shower centroid to the calorimeter 
front face, taken along the direction of its prop¬ 
agation. The algorithm depends on two param¬ 
eters Wq and a, which have been optimized to 
Wq = 4.3 and a = 7 cm by Monte-Carlo simula¬ 
tion and real data from the elastic runs, where 
a 2 mm position resolution (a) at 1.1m and 
4.2 GeV was measured, compatible with the one 
obtained from Monte-Carlo simulations. Posi¬ 
tion resolution when two partially overlapping 
clusters are present, is slightly worse than in the 




FIG. 14. (Color online) Squared missing mass 
associated with the reaction ep —> e'jX for Kin2. 
Total events for Kin2 are represented as inverted 
black triangles, the estimated 7r° contamination is 
represented as green diamonds, the distribution af¬ 
ter the subtraction of accidentals and tt^’s is shown 
as blue open circles. Finally, it is compared with 
the DVCS Monte-Carlo shown as a red solid line. 
In order to remove unnecessary uncertainties due to 
low-missing-mass-squared accidental events, we ap¬ 
ply a cut requiring a missing mass squared higher 
than 0.5 GeV^ for all kinematics. 


case of a single cluster: simulated data show in 
this case a 4 mm spatial resolution. 


C. Event Selection 


The ep —>■ epj events are selected among the 
calorimeter 1-cluster events. A software energy 
threshold of 1.1 GeV was applied to calorimeter 
clusters, slightly above the hardware threshold 
of ~ 1 GeV. Fiducial cuts were used to dis¬ 
card events hitting blocks at the edges of the 
calorimeter. Figure shows the ep —>■ eyV 
missing-mass-squared distribution of the data. 
Accidental coincidences were estimated by ana¬ 
lyzing events in [-11,-5] and [5,11] ns time win¬ 
dows, the same width as the coincidence clus¬ 


tering window but shifted in time (see Fig. 15). 
The use of two intervals to estimate the acciden¬ 
tal sample reduces its statistical uncertainty. 

Neutral pion decays with only one photon 
reaching the calorimeter form an important 
source of background to the DVGS sample. This 
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FIG. 15. (Color online) Time spectrum of blocks 
with E > 300 MeV in the Kin3 setting. It shows 
the 45 ns time window of the waveform analysis. 
The 2 ns CEBAF beam structure is clearly visible. 
The coincidence [-3,3] ns window used for cluster¬ 
ing is shown by the solid line. Dashed lines show 
windows used for the HRS-calorimeter accidental 
subtraction. 


background is subtracted using tt*^ events where 
the two photons are detected in the calorimeter. 
For each detected 7r°, its isotropic decay in its 
center-of-mass frame is simulated Udec = 5000 
times, and the decay photons are projected onto 
the calorimeter acceptance. This simulation al¬ 
lows us to make a statistical subtraction of the 
7 r° background to the DVCS signal, including 
both exclusive and inclusive 7r° events. The sub¬ 
traction is obtained from the simulated decays 
in which one of the photons is emitted close 
to the pion momentum direction. Note that 
this background subtraction scheme could not 
be applied in Kinl as the energy of 7r° decay 
photons is too close to the calorimeter threshold 
to ensure an efficient background subtraction. 


A self-consistency check of the 7r° subtraction 
method was performed using a Monte-Carlo. 
7 r*^’s were generated over the acceptance and 
classified into two categories: the one-photon- 
detected and the two-photon-detected events. 
After applying the 7r°-subtraction method de¬ 
scribed above to the two-photon category, we 
obtained a number of one-photon events and 
compared it to the one-photon-detected cate¬ 
gory. The result is presented on Fig 16 This 
efficiency ratio is close to 1 except in the cor¬ 


ners or close to the edges. Therefore we applied 
a geometrical cut on the cluster in the data and 
the Monte-Carlo simulation, also shown on the 
figure. 



XI (cm) 


FIG. 16. (Golor online) The octagonal cut on the 
front face of the calorimeter applied to all events in 
order to ascertain a high efficiency in the 7r° sub¬ 
traction. The general shape and size of the cut can 
be understood from the size of the front face and 
the width of the shower profile (Fig. [^. 


Figure 17 (top) shows the total distribution 
of events in a Kin3 bin along with the acciden¬ 
tals and 7r° contributions. Accidental events re¬ 
side close to (j) = 0° which corresponds to the 
beamline side of the calorimeter, where higher 
single rates are observed. The contribution of 
7 r° events, however, is larger around (p = 180°. 
This feature remains true for most experimen¬ 
tal bins. The bottom plot of Fig. [17] shows the 
helicity-dependent distribution of events for the 
same bin. The contribution of accidental events 
cancels in this difference of counts, as they are 
essentially helicity-independent. The fact that 
the same feature is observed for tt^’s is not triv¬ 
ial. As it turns out, exclusive 7r° events are 
known to have a small beam-spin asymmetry 
at Jefferson Lab kinematics [M], and the 7r° 
events we subtract may include semi-inclusive 
7 r°’s that have an even smaller asymmetry |55] . 
We have checked that in all our experimental 
bins the contributions of both accidental and 
TT^ events to the difference of counts for oppo¬ 
site helicities are compatible with 0 within sta¬ 
tistical error bars. We have therefore decided to 
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not subtract these contributions in the compu¬ 
tation of the helicity-dependent cross sections. 
In this way, even though we were unable to eval¬ 
uate the unpolarized cross section for Kinl, we 
did succeed in evaluating the helicity-dependent 
cross section. 



FIG. 17. (Color online) Total counts (top) and 
difference of counts for opposite helicities as a func¬ 
tion of 0 (bottom) for the Kind bin xb = 0.37, 
= 2.36 GeV^ and —t = 0.32 GeV^. The solid 
curve histogram in black corresponds to the dis¬ 
tribution of events after all analysis cuts have been 
performed. The estimated contribution correspond¬ 
ing to accidental events is shown as a dashed green 
histogram. The estimated 7r° contribution is repre¬ 
sented as a red dotted histogram. 

After the 7r° subtraction, the only remain¬ 
ing channels (other than exclusive DVCS) are 
kinematically constrained to > {M -l-m 7 r)^. 
However, resolution effects may cause these 
channels to contribute below the cut. This 
contamination was strongly suppressed by the 
tight missing-mass-squared cut and will be eval¬ 
uated in section IVI Al 


HRS 

Track mult. 

R-function 

Cerenkov 

Vertex 

Ntr — 1 

R-value> 0.005 

ADC,„„, > 150 

—6.0 cm< Vz < 7.5 cm 

Calorimeter 

Cluster mult. 

Energy 

Position 
(in cm) 

^clus — 1 

Eclus >1.1 GeV 

— 13 < Xclus < 10 
lyciusl 13 

Ucius < 0.92 X {xcius + 13) + 7.5 
Veins < —0.92 X {Xclus — 4) + 13 
Veins > -0.92 X {Xclns + 13) - 7.5 
Veins > 0.92 X {Xclns 4) 13 

Exclusivity 

All settings 

Kinl-2 

Kin3 

Mi > 0.5 GeVVe* 

Mi < 0.95 GeVVe^ 

Mi < 1.09 GeVVe^ 


TABLE II. Summary of the ep — >■ epy selection 
cuts. 


All selection cuts in this analysis are summa¬ 
rized in Tab. El 


D. Efficiencies and Normalization 

The efficiency of the scintillators that were 
used for the electron trigger was monitored dur¬ 
ing dedicated runs along the experiment. An 
efficiency of 99.95% was measured over the du¬ 
ration of the experiment. The efficiency of the 
Cerenkov counter used to discriminate electrons 
from negative pions was measured to be 99%. 
The purity of the electron sample was estimated 
at 98.8%, further enhanced by the missing- 
mass-squared cut on H(e, e'"f)X. We estimated 
that a maximum of 0.5% of electrons may still 
be misidentified and consider this value as the 
systematic uncertainty on the electron identifi¬ 
cation. 

The dead time associated with the data ac¬ 
quisition is determined by comparing the num- 
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ber of pulses from two clocks running both at 
62.5 MHz: one is always running, the other one 
is veto-ed when the DAQ is busy. The inte¬ 
grated luminosity is corrected for the dead time 
on a run-by-run basis, with an associated sys¬ 
tematic error estimated to be 1% for Hall A [43] . 
The average dead time varied between 14% and 
40% depending on the kinematic setting. 

When multiple tracks were detected in the 
HRS, events were discarded due to the unrelia¬ 
bility of the reconstruction. These events repre¬ 
sent between 7% and 10% of the total statistics, 
depending on the kinematic setting. However, 
most of these multi-track events show a very low 
energy in the pion rejector, indicating that most 
of them contain secondary tracks from show¬ 
ers generated in the exit region of the Q3 mag¬ 
net or pions that trigger the DAQ with 5—rays. 
The number of multi-track events correspond¬ 
ing to good electrons was estimated by requir¬ 
ing a ^1.7 GeV minimum energy deposited in 
the pion rejector. The number of good elec¬ 
tron events with 2 or more tracks in the VDCs 
amounts to only ^2% of the total number of 
events for all kinematics. The 0.5% associated 
systematic error has been evaluated by chang¬ 
ing the energy threshold of the pion rejector. 

Similarly, multi-cluster events in the DVCS 
calorimeter are discarded from the analysis. 
They represent from 1% to 5% of the statistics, 
depending on the kinematic setting. In order 
to apply a correction for this, 2-cluster events 
were thoroughly studied. All selection cuts were 
applied to each of the two photons, and a cor¬ 
rection was computed, based on the number of 
events that remain after the cuts are applied. 
Two-cluster events with an invariant mass be¬ 
tween 100 and 170 MeV/c^ were not included 
in the sample used to calculate this correction 
as they are mostly decay photons from neutral 
pions. In rare cases where both photons fulfilled 
all selection cuts, they contribute to the correc¬ 
tion with a relative weight based on the acci¬ 
dental rate measured in their respective kine- 
matical bin. We attribute a systematic uncer¬ 
tainty to the multi-cluster correction based on 
the number of events with more than 2 clusters 
in the calorimeter, which were not considered 


in our analysis. This number represents 7% of 
the 2-cluster events, and therefore an associated 
systematic uncertainty of less than 0.4% overall. 

Table jTTTj gives a summary of efficiency factors 
applied to experimental yields. 


Source Correction to yield 


Cerenkov 

i.oi 

Multi-track 

1.02 

Multi-cluster 

1.02 

Total 

1.05 


TABLE III. Summary of efficiency factors to be 
applied multiplicatively to experimental yields. The 
multi-cluster correction depends on the kinematic 
setting and the experimental bin, only the average 
value is listed in this table. 


E. Monte-Carlo Simulation 

The experimental setup was implemented in 
a Geant4 Monte-Carlo simulation. The HRS 
geometrical acceptance was modeled by a col¬ 
limator window placed at the entrance of the 
spectrometer. Its acceptance was simulated by 
applying the same R-value cut that is used for 
the experimental data (R-value > 5 mrad). The 
PbF 2 DVCS electromagnetic calorimeter geom¬ 
etry was implemented in detail, including all ac¬ 
tive and passive materials of the experimental 
setup. Only the energy deposit of particles in 
the calorimeter is digitized in our simulation, as 
the generation and tracking of Cerenkov pho¬ 
tons requires unrealistic simulation times and 
proves to be unreliable due to the difficulty to 
define optical surfaces accurately. Detector off¬ 
sets were adjusted following geometrical surveys 
of the experimental equipment. 

Events were generated following a flat distri¬ 
bution in xb, t, (p and (pe- In addition, 
the z-position of the vertex was randomized 
within the full length of the target cell. The 
ranges of and xb are defined by the an¬ 
gular and momentum acceptance of the HRS. 
The hadronic part of the reaction {"f*p —>■ jp) is 
computed in its center-of-mass and final-state 







20 


particles are then boosted to the laboratory 
frame. The generation range in t is kinemati¬ 
cally constrained event-by-event by the values 
of and xb- The angle (f) is then gener¬ 
ated uniformly inside 27r. Finally, all particles 
in the final state are rotated around the beam 
axis by chosen large enough to cover the 
full vertical acceptance of the HRS for all po¬ 
sitions along the length of the target. Each 
event is then weighted by a phase-space fac¬ 
tor AF = Aa;BA(5^A(/)At(a;B, A(/)e/Agen, 
where Ng^n is the total number of generated 
events. 

Because of Bremsstrahlung energy losses 
and resolution effects, the missing-mass-squared 
cut removes a significant fraction of exclusive 
events. This is corrected through the Monte- 
Carlo simulation by applying the same cut in 
the simulated data. However, the experimental 
resolution of the calorimeter and the imperfec¬ 
tions of the calibration procedure have to be 
reproduced by the Monte-Carlo simulation. In 
order to achieve this, the detector is divided into 
49 partially overlapping areas. From the pho¬ 
ton four-momentum in the Monte-Carlo simu¬ 
lation the following smearing transformation is 
applied: 


( 



% 

1 —> gaus{fj,, a) x 

Qy 

Qz 


Qz 

\e) 


\e) 


(23) 


In each area, the parameters /r and <j are fitted 
in order to best match the spectra of the 
simulated and the experimental data in the ex¬ 
clusive region. The final values of /r and a used 
to smear the simulated events are interpolated 
event-by-event according to the impact point of 
the photon in the calorimeter. Figure [T8| shows 
the resulting values of /r and a for Kin3, interpo¬ 
lated across the calorimeter surface, and within 
the fiducial region defined by the octogonal cut 
shown in Fig. 16 The parameter /r corrects im¬ 
perfections in the estimation of the energy in the 
Monte-Carlo simulation compared to the data. 


The parameter a accounts for different resolu¬ 
tions on the different areas of the calorimeter. 
The latter can be due to either different levels of 
background or different quality of the crystals. 



FIG. 18. (Color online) Mean p (top) and standard 
deviation a (bottom) of the gaussian distribution 
used to smear the simulated photon data of Kin3 
viewed on the calorimeter surface ydus vs. x^us 
(beam on the right side). A worse energy resolu¬ 
tion is observed at small angles with respect to the 
electron beam (positive Xcius)- We also notice areas 
of fluctuating resolution corresponding to varying 
quality of the PbF 2 crystals. 


The missing-mass-squared cut to ensure ex¬ 
clusivity is chosen as the value where the Monte- 
Carlo and the data spectra start to differ due 
to contamination by non-exclusive events. This 
leads to two different values of missing-mass- 
squared cut: 0.95 GeV^ for KinI and Kin2, 
1.09 GeV^ for Kin3. A study of the system¬ 
atic uncertainty on the exclusivity is presented 
in section fVl Al 
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F. Cross Section 

To derive differential cross sections from the 
measured data, the solid angle (or acceptance) 
fl of the detection apparatus has to be accu¬ 
rately known. In the expression of the pho¬ 
ton electroproduction cross section, Compton 
Form Factor (CFF) combinations F{Q'^,XB,t) 
appear multiplied by different kinematical fac¬ 
tors T{Q‘^,XB,t,4)), which also vary within the 
bin width. In addition, Bethe-Heitler itself is 
a rapidly varying cross section, especially as a 
function of 0, xb and t. Since all the kinematic 
dependences besides the intrinsic CFF ones are 
known, we decided to use a method which di¬ 
rectly extracts the CFF from data by disentan¬ 
gling all effects in a combined data-Monte-Carlo 
fit. This method has the additional advantage 
of automatically handling bin migration effects 
that may occur. The extraction method is for¬ 
mally described in the following. 

Let 


( Eb\ 

Xb 

t 


\ (jbe J 


(24) 


represent the kinematic variable vector at the 
vertex in the simulation. The incident electron 
energy Eb is included in order to treat the ra¬ 
diative tail. Let 


/ Eb\ 


Xb 



0 

\^e J ^ 


(25) 


represent the reconstructed event variables. In 
the Monte-Carlo simulation, we define the map¬ 
ping 


(26) 


as the conditional probability distribution to 
observe an event at the kinematic point Xe 
starting from vertex point x.^. The experimen¬ 
tal acceptances, intrinsic detector efficiencies 
and resolutions, and real-radiative effects are 
included in K{xe\'x.y). This conditional proba¬ 
bility, which we compute using the Monte-Carlo 
simulation, takes into account the potential bin 
migration due to detector resolutions and radia¬ 
tive effects. The binning vector 


i.= ’«■ (27) 

V 4 / e 

labels a set of bins in the corresponding event 
kinematics, after integration over (j)f, since un¬ 
polarized target observables depend on only one 
azimuthal angle (j). The binning vector 


j 


V — 


^ jxB \ 

V jt J 


(28) 


labels a similar set of bins in the vertex vari¬ 
ables. The helicity-dependent and helicity- 
independent cross sections can be written as 
a sum of several harmonic contributions as de¬ 
scribed by Eq. (|4]|^: 

a(x,) = ^r^(x,)xA, (29) 

A 


where F''^(x^) represent some kinematical fac¬ 
tors and XA are some combinations of CFFs 
that are unknown and that parametrize the 
DVCS cross section. Notice that, as shown by 
Eq. (281, the variable (j) is not binned at the 
vertex. This is because the full (/(-dependence 
of the cross section is known and contained in 
the kinematical factors F^(x.u). Thus, the un¬ 
knowns XA are independent of (j). The fact that 
the total number of bins in the reconstructed 
event variables is significantly higher than the 
number of bins in the vertex variables is pre¬ 
cisely what makes the fit described below pos¬ 
sible. 


iL(Xe|Xi,) 
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The number of counts per bin at the vertex is 

iv(j„) = c f Yl = 

Jx„eBin(j„) ^ 

[ r^(x„)rfx„, (30) 

Jx„eBin(j„) 

where C is the integrated luminosity. In the 
experimental bin ig, the yield is 


m^e) 


= j dXe^ A^(i«)ii:(Xe|x„) 

XeeBin(ie) J” 


j dXe J dx^r^{Xy)K{lCe\Xy) 

j” XeeBin(ie) x„eBin(j„) 

(31) 


We define a bin mapping function: 


J dXedXy K{Xe\Xy)T^{Xy) . 
XeeBin(ie) x„eBin(j„) 

(32) 

This function is basically the solid angle 
weighted by the kinematic factors r^(x^,), 
where the effects of bin migration are taken into 
account through the function K{xe\xy). Gener¬ 
ally, the number of counts per bin can thus be 
written as: 

N^^{ie) = ■ (33) 

j..A 


Note the summation over all jy. All bins at the 
vertex might contribute to a given experimental 
bin ie, with a certain probability or weight given 
by the function ATA ^ , computed in the simula¬ 
tion. We construct a which we minimize to 
extract the Xj^: 

, ^ [jVExp(i^)_jVMC(i^)]2 

A [5Exp(i^)]2 ’ 


where i5®’'P(ie) are the experimental statistical 
uncertainties in each bin. 

The coefficients Xj^ are defined as the values 


of Xj^ that minimize x^- 


0 = - 


1 dx'^ 


2 5XA 


0 = E «jt j' EJ ^ j- A • (35) 

It. A' 

The linear system is defined by: 


A.A' 

«j.jt=E^ Ik 


/3a = E^ 


[kExp(i^)]2 ’ 

Af^^’^P(ie) ATA 


[<5Exp(i^)]2 


The fit parameters are: 
^A 


4 = X ^ 

Ik A' 


(36) 

(37) 

(38) 


The covariance matrix of the fitted parameters 
is: 


^A,A' 

' Jt 


^ A' 


(39) 


Finally, the cross-section values (and associ¬ 
ated error bars) at the point Xi^ are obtained 
as: 


dV(5:iJ _ N^^P{iy,iy) 

dxsdQ'^dtdcj) dxsdQ'^dtdcj) ’ 

(40) 
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where 


(a;ie) 


dxsdQ'^dtdcj) 


= (41) 


is defined by the fit parameters of the bin 
which has the same bin limits in xb, and 
t as the experimental bin ig. The number of 
counts -/V“‘"(ie,jt,) and j„) corrected 

from bin migration are given by: 


N 


MC, 




(42) 


-^E E <jt^E(43) 

A 

Note that we study the harmonic coefficients 
of the cross section which are the sum of (c, s)^ 
and themselves involving several 

combinations of twist-2 and twist-3 CFFs. The 
(()-dependence is therefore not enough to sepa¬ 
rate all linear/bilinear combinations of CFFs. 
The 0-dependence of the cross section can thus 
be properly described by different choices of free 
parameters. 

In this analysis, we chose to parametrize the 
DVCS helicity-independent cross section by the 
three following combinations of effective CFFs: 
CDVCS^jr^j:*'^ (Ec^, 3fie[C^(J')] (Eq. and 
3?e[C^(J^e//)] (Eq. 111. The helicity-dependent 
cross section is fitt^using the 5m[C^(J^)] and 
9m[C^(Je//)]. Three reasons have led to this 
choice: 

• The contributions to the cross section as¬ 
sociated to each of these parameters have 
a distinct 0-dependence, minimizing the 
correlations among them, 

• We keep the dominant twist DVCS^ con¬ 
tribution, 

• Higher twist contributions are kinemati¬ 
cally suppressed. 

While this is the most physical choice of pa¬ 
rameters, any other choice that provides a good 
fit {x^ /dof ~ 1) to the 0-dependence of the 


number of counts is an equally valid choice as 
far as the cross-section extraction is concerned. 
The fitted parameters, though, would have a 
less straightforward physics interpretation in 
that case. We have tested the stability of our 
cross-section results against a different choice 
of free parameters and results are discussed in 
section IVI Al 


G. Global Normalization 


In the previous section, we defined C as the 
integrated luminosity. It is computed from the 
average total charge Q recorded by the BCMs 
as: 


J dt e Ah ' 


(44) 


where e = 1.602 • 10“^® C is the electron charge. 
Ah = 1.0079 g/mol is the atomic mass of H, and 
Na = 6.022 • 10^^ mol“^ is Avogadro’s number. 
The LH 2 target length was I = 15 cm and was 
operated at 19 K and a pressure of 25psi, which 


shows the integrated luminosity (corrected by 
the acquisition deadtime) recorded for each of 
the kinematic settings. 


gives a density of p = 0.07229g/cm^. Table IV 


Kin 

Q+(C) 

Q-iC) 

0(C) 

Qasy (10-3) 

£ (fb-3) 

1 

0.3732 

0.3733 

0.7464 

-0.1 

3059 

2 

0.4057 

0.4064 

0.8121 

-0.7 

3328 

3 

0.6913 

0.6937 

1.385 

-2.4 

5676 


TABLE IV. Helicity-correlated charge (Q±), total 
charge {Q = Q+ -I- Q-) and charge asymmetry 
Qaay = {Q+ — Q -) jQ for the three kinematics set¬ 
tings. The last column shows the integrated lumi¬ 
nosity including events for which the helicity bit is 
undefined. In all cases, the charge and the luminos¬ 
ity have been corrected for the deadtime. 


V. RADIATIVE CORRECTIONS 

The diagrams of Fig. [^only include the low¬ 
est order QED amplitude for the DVCS process. 
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The experimental cross section necessarily in¬ 
cludes higher-order QED processes. Radiative 
corrections to the ep —> epy reaction have been 
studied in several papers [5SH5U] . In this anal¬ 
ysis we follow the approach of Hg. 

In this analysis we fit a model cross section to 
the experimental yield, bin by bin. We separate 
the radiative corrections into terms that are de¬ 
pendent on the missing-mass-squared M|- cut 
we impose on the H(e, e'j)X spectra, and terms 
that are independent of this cut. The external 
radiative effects on the incident electron, and 
internal real radiative effects at the vertex are 
treated in the equivalent radiator approxima¬ 
tion isoiini]. Pre-scattering radiation is mod¬ 
eled by generating an event-by-event energy loss 
Aifin of the incident electron (Ei,) following a 
distribution {b ~ 4/3): 


liniEb, XEin.tin) = 


bt in + <^s/2 

AEin 


^Ein 

Eb 


with 


TT 


mi 


bti„-|-<5s/2 

(45) 

(46) 


where tin is the event-by-event target thickness 
(in radiation lengths) traversed by the electron 
before the scattering vertex. The Schwinger 
term 5s models the internal pre-scattering ra¬ 
diation. The scattered energy at the vertex is 
E'^ = Eb — AEin — I{2MpXB)- Internal post¬ 
scattering radiation is modeled by a similar dis¬ 
tribution in the post-scattering radiated energy 
AEont'- 


Ss/2 

Aifout 


AElout 

. E', 


i5s/2 


(47) 


These radiative effects are treated within 
the peaking approximation. External post¬ 
scattering radiation by the scattered electron is 
evaluated with the Monte-Carlo simulation by 
transporting the electron to the entrance of the 
spectrometer. Kinematic shifts {e.g. in either 
the norm or the direction of q) from external 
and internal radiations are fully included in the 
simulation and thereby unfolded from the ex¬ 
tracted cross sections. 


In addition to these radiative effects incor¬ 
porated into our Monte-Carlo, we correct the 
data for internal virtual radiation as well as the 
cut-off independent effect of unresolvable soft 
real radiation, given by Eqs. 58-62 of [Sg. The 
virtual corrections to the VCS amplitude are 
model independent, in the sense that they do 
not depend on the dynamics of the 7*p —>■ yp 
process. These corrections (vacuum polariza¬ 
tion and vertex renormalization) are essentially 
equivalent to the corrections to elastic ep scat¬ 
tering, with suitable adjustment to the kinemat¬ 
ics. On the other hand, the vacuum polariza¬ 
tion and vertex corrections to the BH ampli¬ 
tude differ by several percent relative to the 
VCS corrections, and the BH amplitude also 
has self-energy corrections to the virtual elec¬ 
tron propagators. We calculate separately the 
radiative corrections to the helicity-independent 
and helicity-dependent cross sections based on 
a code derived from |56] which includes the 
leading-twist DVCS amplitude with a fully fac¬ 
torized GPD ansatz [Sg. The correction 
factors vary by less than 0.5% over </) and by 
« 1% over the [xs, acceptance of each kine¬ 
matic setting. We assign a 2% systematic er¬ 
ror to the combined real- and virtual-radiative 
corrections. This is based on the variation of 
the correction over the acceptance, ambiguities 
over whether or not to exponentiate the cor¬ 
rection, and the model-dependence of the rel¬ 
ative contributions of the |T^^| , interference 

and terms in the unpolarized cross 

sections. Over our five kinematic settings, the 
average corrections varied by less than 0.5%. 
Since this is less than the uncertainty of the 
correction, we apply the following global cor¬ 
rections to all cross-section bins: 

^4^Born ^ (0.948 ± 0.02) 

A4^Born ^ (q gyg ^ q q 2) A^cr^’^p. (48) 

VI. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are divided into un¬ 
correlated (or point-to-point) and correlated 
(or normalization) uncertainties. The largest 
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source of uncorrelated error in this experiment 
was associated to the missing-mass-squared cut. 
The correlated uncertainties have been de¬ 
scribed before and a summary table is shown 
in this section. 


A. Missing-mass-squared Cut 


Two systematic effects are associated with 
the missing-mass-squared cut. The first comes 
from semi-inclusive events contaminating our 
sample. These events have larger missing-mass- 
squared values induced by extra missing parti¬ 
cles. Indeed, even if the cut is supposed to keep 
this contamination minimal, a small fraction 
of such events may remain below the missing- 
mass-squared cut. In order to evaluate an up¬ 
per value for this systematic error, we exam¬ 
ined the ratio of the integrals of the experi¬ 
mental and Monte-Carlo missing-mass-squared 
spectra. As seen in Fig. this ratio increases 
significantly with the missing-mass-squared cut, 
which is expected since the Monte-Carlo only 
contains exclusive events. By varying the cut 
from the nominal value 0.95 GeV^/c^ up to 
1 GeV^/c"^, the observed contamination remains 
smaller than 1%, which we took as the system¬ 
atic uncertainty on the cross section. 

The second effect induced by the missing- 
mass-squared cut arises from a mismatch on the 
position and shape of the missing-mass-squared 
peaks between data and Monte-Carlo. This is 
due to our limited ability to reproduce perfectly 
the response of our calorimeter. This mismatch 
increases as the missing-mass-squared cut de¬ 
creases and is maximal around the maximum of 
the distribution. We estimate the correspond¬ 
ing error by looking at the variation of the cross 
section between the nominal cut and a lower cut 
value. This lower bound is chosen such that the 
loss of statistics is 15%, ensuring that the ob¬ 
served variations are not statistical in nature. 
The systematic error is evaluated for each (t, (p) 
bins of each kinematic setting and may reach up 
to a few percent. These point -to-point uncer¬ 
tainties are included in the data tables IVlHfXVI 



FIG. 19. (Color online) Top: Variation of the 
ep —>■ ep 7 cross section for Kin2, —t — 0.17 GeV^, 
as a function of the missing-mass-squared cut, for 
<j) = 0° (upper blue points) and p = 180° (lower 
black points). The dotted vertical line corresponds 
to the nominal cut. The systematic errors are eval¬ 
uated bin by bin in p and t for each kinematic set¬ 
ting by studying the variation of the cross section 
between the nominal and the lower missing-mass- 
squared cut (dashed line). The insert represents the 
same cuts on the missing-mass plot. Bottom: Ra¬ 
tio of the integrals of the experimental and Monte- 
Carlo missing-mass spectra, as a function of the 
missing-mass-squared cut. By varying the cut up 
to 1 GeV^/c'^ , represented by the dotted-dash line, 
the observed contamination remains smaller than 
1 % (green band). 


B. Cross-section Parametrization 

As mentioned in section EVFl the cross- 
section results should be independent of the 
choice of parametrization in the extraction 
method. To evaluate the impact of this choice, 
we used a different parameter set by replacing 
the squared DVCS amplitude term by the inter¬ 
ference term 3?e , which yields an equally 

good fit to the data. A difference in the cross- 
section value of up to 1% appears locally de¬ 
pending on the kinematic bin, as shown in 
Fig[^ As a consequence we estimated the sys- 
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tematic error from the parameter choice to be 

1 %. 



FIG. 20. Difference in % between the cross sec¬ 
tion extracted with the squared DVCS amplitude 
term and with the 5Re C^’'^ term for xb = 0.37, 
= 2.36 GeV^ and -t = 0.33 GeV^ The 0- 
profile of the difference is a consequence of the small 
COS0 and cos2(/) dependences of the 5Re kine¬ 
matic coefficient. Both extractions give almost the 
same reduced /dof=0.94 (nominal) and 0.93 (al¬ 
ternate) for the entire Kin2 setting. 


C. Correlated Uncertainties 

Table presents the systematic uncertain¬ 
ties on the cross section stemming from nor¬ 
malization effects, which are considered 100% 
correlated bin-by-bin. Note that the helicity- 
dependent cross sections have an additional un¬ 
certainty coming from the beam polarization 
measurement. The determination of these un¬ 
certainties are discussed in the associated sec¬ 
tion listed in the table. 


VII. RESULTS 

The cross-section extraction procedure de¬ 
scribed in section HYl was applied to all 
data sets, for both the unpolarized and the 
helicity-dependent cases. In addition to the 
Q^-dependence of the helicity-dependent cross 


Systematic uncertainty Value Section 


HRS acceptance cut 

1% 

Electron ID 

0.5% 

HRS multitrack 

0.5% 

Multi-cluster 

0.4% 

Corrected luminosity 

1% 

Fit parameters 

1% 

Radiative corrections 

2% 

Beam polarization 

2% 


IV A 
IV D 
IV D 
IV D 
IV D 
VIB 

[v 

III A 3 


Total (helicity-independent) 2.8% 
Total (helicity-dependent) 3.4% 


TABLE V. Normalization systematic uncertainties 
in the extracted photon electroproduction cross sec¬ 
tions. The systematic error coming from the fit pa¬ 
rameter choice is not a normalization error per se, 
but we consider that 1% is an upper limit for this 
error on all kinematic bins. The helicity-dependent 
cross sections have an extra uncertainty stemming 
from the beam polarization measurement. The last 
column gives the section in which each systematic 
effect is discussed. 


sections, we were able to measure the Q^- 
dependence of the unpolarized cross section at 
two values of Q^=1.9 and 2.3 GeV^. The 
xs-dependence of helicity-dependent and - 
independent cross sections were studied using 
the KinX2 and KinX3 settings. Note that an 
extra bin in t was analyzed compared to our 
previous publication [29] for all {xb,Q^) set¬ 
tings. 

An example of the cross-section extraction 
is presented in Fig. ^ for xs = 0.37, = 

2.36 GeV^ and —t = 0.32 GeV^, along with 
the different contributions resulting from the fit, 
which gave an overall /dof of 1.1. For the 
unpolarized cross section, one observes a sig¬ 
nificant contribution from the term associated 
with \rDvcsy ’, in addition to a large contribu¬ 
tion from the interference term. Note that the 
YfDVCS I contribution is (()-independent in con¬ 
trast to the BH and interference contributions. 
Indeed, it does not contain the 'Pi{(j))V2{4’) elec¬ 
tron propagators as shown in Eq. [^ The pre¬ 
cision of the data is such that other contribu¬ 
tions than the Bethe-Heitler are obviously nec- 
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<t> [deg] 


FIG. 21. (Color online) Unpolarized (top) and 
helicity-dependent (bottom) cross-section extrac¬ 
tion for the Kind bin —t = 0.32 GeV^. The error 
bars on the data points are statistical only. The 
shaded areas represent the statistical uncertainty 
for each contribution. 


essary to explain the observed cross section. 
The helicity-dependent cross section is domi¬ 
nated by the twist-2 interference term, as no¬ 
ticed before in this experiment [21] and else¬ 
where [2H1ES] ■ These conclusions extend to all 
bins in our an alysi s, w hose results are shown in 
Tab. VI lists the x^/do/ resulting 


VIID 


section 

from the extraction method for all kinematics 
settings. 


A. Scan in 

The combinations of effective CFFs which 
have been extracted from the fitting procedure 
for Kinl-3 using the formalism developed in |40j 
are shown integrated over t in Fig.[^ With the 


Settings 

xloi/dof 

)Cunp/dof 

Kinl 

0.88 

- 

Kin2 

1.00 

1.16 

KinX2 

0.96 

0.82 

Kin3 

1.15 

0.99 

KinX3 

1.08 

1.28 


TABLE VI. /dof resulting from the extraction 
method for all kinematics settings. The subscript 
”pol” stands for polarized cross sections, ”unp” for 
unpolarized cross sections. 


choice of parameters used to describe the kine- 
matical dependen ce of t he cross sections (as ex¬ 
plained in section IV FI, the contribution asso¬ 
ciated with the jg large for the 

unpolarized case. The twist-2 interference term 
is significant and the contribution of the twist- 
3 interference term is often found to be small, 
with large systematic errors. For the polar¬ 
ized case, the twist-2 interference term is dom¬ 
inant, the twist-3 contribution is small, again 
with large systematic errors. 

Overall, the extracted parameters show 
no -dependence for either the helicity- 
dependent or the helicity-independent cases 
over our Q^-range. Note that the logarithmic 
Q^-evolution can safely be neglected within this 
lever arm at this xb- 

The full set of results for settings Kinl-3 are 


presented in Fig. 24 -28 in section VIID 


B. Scan in xb 

The results from KinX2 and KinX3 showing 
the X B-dependence of the cross sections are pre¬ 
sented in Fig. |29| - |32| in section |VII D[ KinX3 
has a limited acceptance close to 0°, which 
increases the correlation between the different 
fit parameters describing the azimuthal depen¬ 
dence of the cross section. Indeed, the separa¬ 
tion of the real part of the twist-2 interference 
and | 7 "-DVCS| contributions in the fit is par¬ 
ticularly sensitive to the relative value of the 
cross section measured around both (j) = 0 and 
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1.4 1.6 1.8 2 2.2 2.4 1.4 1.6 1.8 2 2.2 2.4 

[GeV^] [Qgv^] 

FIG. 22. Combinations of effective CFFs extracted from the fitting procedure described in section [IV F| 
using the formalism developed in [40], integrated over t and plotted as a function of Q^. The top three 
plots show the effective CFFs resulting from the unpolarized cross section fit (Kin2 and Kin3), whereas the 
bottom plots show the effective CFFs resulting from the helicity-dependent cross section fit (Kinl-3). The 
shaded areas represent systematic errors. 


180°. These difficulties have basically no im¬ 
pact on the determination of the cross sections 
themselves. The measured xs-dependence will 
set interesting constraints on GPD models and 
parametrizations, especially thanks to the rela¬ 
tively high accuracy of our data. 

C. Comparison with GPD Models 

In Fig. |23[ we compare our results with vari¬ 
ous models and previous fits to data. We have 
chosen to use two different kinds of double¬ 
distribution GPD models, namely the VGG |M| 
and KMS12 [551 models. Note that in contrast 
to VGG, the KMS12 model was tuned using 
vector meson data at low to very-low x b , and is 
not considered adapted yet to the valence quark 


region. In any case, one observes that both 
models overshoot the helicity-dependent cross 
section data in this Kin2 bin, whereas VGG is 
more adequate for the unpolarized data. 

In addition, we have compared our data with 
the KMlOa model [66|, which fits some of its 
parameters to all DVGS data available world¬ 
wide except for the previously published results 
from a subset of the present experiment. The 
consequence is that no absolute DVGS cross- 
section data in the valence region were used for 
this fit. The KMlOa model is clearly very close 
to the helicity-dependent data, which is not a 
surprise considering that the GLAS asymmetry 
data in the same kinematic region were used to 
constrain this model. However, this same model 
significantly underestimates the DVGS unpolar¬ 
ized cross section around (j) = 180°. 
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Recently, kinematic twist-4 target-mass and 
finite-t corrections (TMC) have been calculated 
for DVCS on the proton and estimated for the 
KMS12 model [571115] (shown in Fig. 23). Since 
this model is not adapted to the valence quark 
region, we have extracted the correction fac¬ 
tor and applied it to the KMlOa parametriza- 
tion [55]. This allows us to gauge the effect 
of such corrections in the most realistic model 
available to us. It is striking that the lack of 
strength observed at (j) = 180° for the KMlOa 
model is largely compensated by the TMC, giv¬ 
ing a surprisingly good agreement between this 
modified KMlOa model and our data. 

An update of the KMS12 model, taking into 
account the DVCS data in the valence region, 
would allow for a much stronger statement 
about the necessity of target-mass and finite- 
t corrections at these moderate Q^. At any 
rate, we emphasize that the high accuracy of 
the present data is crucial to disentangle the dif¬ 
ferent contributions at play in this critical area 
around 180°. There is no doubt that the addi¬ 
tion of our new data set to the KM fit will be 
most interesting, especially in the light of these 
new higher-twist calculations. 

All the features we have described remain 
true for most of our data bins, which are shown 
in section |VII D| It is interesting to note that 
for the highest bins in t, especially for Kin2 
and KinX2 (Figure 24 and 29), the TMC to 
the unpolarized cross section is of the same or¬ 
der as the cross section itself around (j) = 180°. 
This corresponds to values of ~ 0.15 

or larger. It is not unreasonable to expect that 
higher-order corrections in start to be 

important at these values, and may compen¬ 
sate the peculiar behavior of the TMC around 
(f) = 180°, which is not visible in data. Efforts 


2\k 


se- 


to achieve a resummation of the {—tjQ' 
ries to all orders are currently undertaken [70]. 


D. Results for All Kinematics 

In the following we present the unpolarized 
cross sections for Kin2, Kin3 as well as KinX2 



FIG. 23. (Color online) Unpolarized (top) and 
helicity-dependent (bottom) cross sections for the 
Kin2 bin —t = 0.23 GeV^. The light blue area rep¬ 
resents the point-to-point systematic uncertainties 
added linearly to the normalization error. The pre¬ 
dictions from the distribution-based models KMS12 
and VGG are shown as the dashed green and solid 
red curves, respectively. The KMlOa fit is rep¬ 
resented as the solid blue line. The target-mass 
and finite-t corrections are included in the KMS12 
model and shown as the dotted-dash curve. The 
correction is then applied to the KMlOa model 
shown as the dotted blue line. 


and KinX3 in Fig. and [^ respec¬ 

tively, for a total of 468 experimental bins in 
The cross-section differences 
for opposite beam helicities are presented for 
Kinl-3, KinX2 and KinX3 in Fig. [MH28| and 
[3l] [3^ for a total of 588 experimental bins 
in (xb, —t, </>). All results are compared to 
only two models for clarity: the KMlOa model 
and its modihed version, including the TMC ef¬ 
fects as described in section lVII Cl All the cross- 
section data are also listed in Tables IVIIHXVl 
along with their statistical and point-to-point 
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systematic uncertainties. The correlated sys¬ 
tematic errors are summarized in Tab. El 

VIII. SUMMARY AND CONCLUSION 

We report final results for EOO-110, the first 
dedicated Deeply Virtual Compton Scatter¬ 
ing experiment, which ran in Hall A of Jef¬ 
ferson Lab. Using new developments in the 
parametrization of the DVCS reaction, we ex¬ 
tracted cross sections including for the first 
time an evaluation of the DVCS squared am¬ 
plitude. We showed results for the unpolar¬ 
ized DVCS cross section at two different Q^- 
values and two different xs-values, thanks to a 
new analysis which allowed for a reliable eval¬ 
uation of the 7r° background in all these kine¬ 
matics. The effective Compton Form Factors 
used to described the kinematical dependence of 
the helicity-dependent and helicity-independent 
cross section show no -dependence, compat¬ 
ible with the dominance of the leading-twist di¬ 
agram in this region of moderate and high 
xb- Our results were compared with various 
models based on the Generalized Parton Dis¬ 
tributions framework. A relative good agree¬ 
ment was found with the KMlOa parametriza¬ 
tion. However, this model does not fully match 
the behavior of the unpolarized cross section for 
(j) ~ 180°. We showed that adding an empirical 
estimate of the target-mass and finite-f correc¬ 
tions to the KMlOa model improved the agree¬ 
ment with our data significantly, which may 
hint at the necessity to include such effects in 
the analysis of moderate-Q^ data, highly rele¬ 
vant for current and future Jefferson Lab ex¬ 
periments. At any rate, the accuracy of the 
unpolarized cross-section data around (j) = 180° 


seems absolutely critical to disentangle all con¬ 
tributions of the cross section. 

The significant deviation of the DVCS cross 
section observed in this experiment with re¬ 
spect to the Bethe-Heitler contribution moti¬ 
vated the subsequent experiment E07-007 [3T], 
currently under analysis. Its goal is to investi¬ 
gate the nature of this deviation by using the 
beam-energy dependence of the different terms 
of the cross section. Indeed, the BH-DVCS in¬ 
terference contribution has a ~ dependence 
whereas the DVCS^ varies as ^ In a way 
similar to a Rosenbluth separation, by measur¬ 
ing the DVCS cross section at exactly the same 
kinematics but different beam energies, one will 
be able to tell if this deviation is mainly due to 
the DVCS^, the BH-DVCS interference terms 
or higher twist terms. This high accuracy mea¬ 
surement of the cross section at two beam ener¬ 
gies will set stringent constraints on GPD mod¬ 
els. 
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( deg ) 

xg = 0.343 

Q2 = 1.820 GeV ^ 
t = - 0.172 GeV ^ 

Xg = 0.368 
= 1.933 GeV ^ 
t = - 0.232 GeV ^ 

Xg = 0.375 
= 1.964 GeV ^ 
t = - 0.278 GeV ^ 

xq = 0.379 
= 1.986 GeV ^ 
t = - 0.323 GeV ^ 

Xg = 0.381 
= 1.999 GeV ^ 
t = - 0.371 GeV ^ 

7.5 

111.6 ± 4.1 

- 1 - 2.3 

- 0.0 

76.9 ± 4.2 

- 1 - 0.0 

- 1.1 

75.2 ± 5.6 

- 1 - 2.2 

- 1.7 

66.2 ± 7.7 

- 1 - 3.2 

- 0.9 

70.6 ± 11.8 

+ 13.1 

- 0.0 

22.5 

117.6 ± 3.9 

- 1 - 0.7 

- 0.3 

68.6 ± 3.7 

- 1 - 0.0 

- 2.8 

59.7 ± 4.5 

-f 1.4 

- 0.6 

64.1 ± 6.3 

- 1 - 0.9 

- 2.0 

66.3 ± 8.7 

- 1 - 5.2 

- 0.4 

37.5 

99.8 ± 3.5 

- 1 - 0.3 

- 0.5 

61.7 ± 3.2 

- 1 - 0.0 

- 2.6 

57.9 ± 3.8 

- 1 - 0.9 

- 1.5 

43.9 ± 4.5 

-t 1.8 

- 0.0 

45.1 ± 5.7 

- 1 - 0.0 

- 3.3 

52.5 

94.1 ± 3.3 

- 1 - 0.0 

- 3.4 

66.9 ± 3.0 

- 1 - 0.0 

- 1.2 

50.4 ± 3.2 

- 1 - 0.0 

- 4.2 

42.5 ± 3.8 

- 1 - 0.0 

- 2.8 

39.0 ± 4.5 

- 1 - 0.0 

- 5.8 

67.5 

88.2 ± 3.0 

- 1 - 0.0 

- 2.1 

55.5 ± 2.7 

- 1 - 0.0 

- 2.0 

40.0 ± 2.8 

- 1 - 0.2 

- 1.1 

35.3 ± 3.1 

- 1 - 0.0 

- 2.9 

41.0 ± 3.7 

- 1 - 0.0 

- 2.4 

82.5 

78.7 ± 2.8 

- 1 - 0.4 

- 0.4 

46.1 ± 2.4 

-t 0.1 

- 0.9 

36.1 ± 2.5 

-t 1.3 

- 0.0 

29.4 ± 2.8 

-f 1.2 

- 0.4 

26.7 ± 2.9 

- 1 - 0.0 

- 3.3 

97.5 

67.2 ± 2.6 

- 1 - 3.6 

- 0.0 

40.1 ± 2.1 

-t 0.1 

- 1.4 

35.2 ± 2.2 

- 1 - 0.6 

- 1.8 

27.4 ± 2.4 

- 1 - 0.0 

- 2.0 

29.2 ± 2.8 

- 1 - 0.0 

- 3.3 

112.5 

60.8 ± 2.4 

- 1 - 0.5 

- 1.3 

37.7 ± 2.0 

- 1 - 0.0 

- 3.1 

33.0 ± 2.1 

- 1 - 0.0 

- 1.4 

24.0 ± 2.1 

- 1 - 0.0 

- 1.3 

25.9 ± 2.6 

- 1 - 0.0 

- 1.0 

127.5 

57.5 ± 2.3 

- 1 - 0.0 

- 2.3 

34.0 ± 1.8 

- 1 - 0.0 

- 2.2 

29.3 ± 1.8 

-t 0.1 

- 3.0 

25.5 ± 2.0 

- 1 - 0.0 

- 2.2 

21.4 ± 2.4 

-t 0.1 

- 1.4 

142.5 

50.1 ± 2.1 

- 1 - 0.2 

- 0.8 

33.9 ± 1.7 

- 1 - 0.0 

- 0.9 

28.0 ± 1.7 

- 1 - 0.0 

- 1.4 

22.4 ± 1.9 

- 1 - 0.5 

- 1.3 

16.1 ± 2.4 

- 1 - 0.0 

- 1.4 

157.5 

49.9 ± 2.1 

- 1 - 0.0 

- 0.9 

34.6 ± 1.7 

- 1 - 0.2 

- 0.5 

26.3 ± 1.7 

- 1 - 0.0 

- 2.8 

21.3 ± 2.1 

- 1 - 0.0 

- 2.6 

13.5 ± 2.7 

- 1 - 2.7 

- 0.0 

172.5 

48.1 ± 2.0 

- 1 - 0.0 

- 2.0 

30.7 ± 1.6 

- 1 - 0.9 

- 0.3 

26.8 ± 1.7 

- 1 - 0.0 

- 0.7 

17.4 ± 2.0 

-t 1.3 

- 0.0 

17.5 ± 3.2 

- 1 - 0.0 

- 1.7 

187.5 

48.3 ± 2.0 

- 1 - 0.0 

- 1.5 

31.5 ± 1.6 

- 1 - 0.0 

- 1.6 

26.6 ± 1.7 

- 1 - 0.2 

- 0.6 

21.4 ± 2.1 

- 1 - 0.7 

- 0.3 

20.5 ± 3.2 

- 1 - 0.0 

- 1.4 

202.5 

53.2 ± 2.1 

- 1 - 0.0 

- 1.6 

35.7 ± 1.7 

-f 1.2 

- 0.0 

24.7 ± 1.6 

- 1 - 0.0 

- 1.5 

22.5 ± 2.1 

- 1 - 0.6 

- 0.3 

21.7 ± 2.9 

- 1 - 0.0 

- 1.8 

217.5 

52.7 ± 2.2 

- 1 - 0.0 

- 0.7 

33.3 ± 1.7 

- 1 - 0.6 

- 0.0 

28.4 ± 1.8 

- 1 - 0.0 

- 4.2 

25.1 ± 2.1 

- 1 - 0.0 

- 1.7 

26.5 ± 2.6 

- 1 - 0.0 

- 1.4 

232.5 

55.2 ± 2.2 

- 1 - 0.0 

- 3.4 

33.7 ± 1.8 

- 1 - 0.2 

- 1.1 

28.4 ± 1.9 

- 1 - 0.0 

- 1.5 

26.0 ± 2.1 

- 1 - 0.0 

- 0.4 

20.1 ± 2.2 

- 1 - 0.0 

- 2.4 

247.5 

58.5 ± 2.3 

- 1 - 0.2 

- 0.9 

37.7 ± 2.0 

- 1 - 0.0 

- 1.3 

37.3 ± 2.2 

- 1 - 0.0 

- 2.5 

30.5 ± 2.4 

-t 0.1 

- 3.1 

25.7 ± 2.6 

- 1 - 0.0 

- 1.8 

262.5 

63.1 ± 2.4 

- 1 - 0.0 

- 2.8 

40.4 ± 2.1 

- 1 - 0.0 

- 1.7 

37.8 ± 2.3 

- 1 - 0.0 

- 2.5 

28.2 ± 2.5 

-t 1.3 

- 0.3 

21.7 ± 2.6 

- 1 - 0.6 

- 1.6 

277.5 

74.2 ± 2.7 

- 1 - 0.8 

- 0.7 

42.9 ± 2.3 

- 1 - 0.2 

- 0.7 

40.2 ± 2.4 

- 1 - 0.0 

- 2.6 

32.6 ± 2.7 

- 1 - 0.0 

- 1.7 

26.2 ± 2.9 

- 1 - 0.5 

- 0.3 

292.5 

81.8 ± 2.9 

- 1 - 0.3 

- 2.1 

52.8 ± 2.6 

- 1 - 0.6 

- 1.2 

42.6 ± 2.7 

-f 1.1 

- 0.1 

36.4 ± 3.1 

- 1 - 3.0 

- 0.0 

29.7 ± 3.5 

- 1 - 0.4 

- 0.7 

307.5 

96.7 ± 3.2 

- 1 - 0.0 

- 1.8 

59.3 ± 2.8 

- 1 - 0.0 

- 2.6 

46.8 ± 3.1 

- 1 - 0.4 

- 0.4 

47.4 ± 3.8 

- 1 - 0.0 

- 3.2 

33.8 ± 4.1 

- 1 - 0.0 

- 4.1 

322.5 

105.6 ± 3.5 

- 1 - 0.0 

- 1.2 

62.0 ± 3.2 

- 1 - 0.0 

- 1.9 

48.0 ± 3.5 

- 1 - 0.0 

- 2.0 

46.0 ± 4.7 

- 1 - 0.7 

- 1.8 

39.6 ± 5.7 

- 1 - 0.0 

- 5.5 

337.5 

108.6 ± 3.8 

-t 1.5 

- 1.5 

73.3 ± 3.8 

-f 4.1 

- 0.0 

60.3 ± 4.5 

-t 1.6 

- 0.5 

52.6 ± 6.3 

-f 1.4 

- 0.8 

44.9 ± 8.2 

- 1 - 6.9 

- 0.0 

352.5 

113.6 ± 4.0 

- 1 - 2.9 

- 0.0 

76.8 ± 4.2 

- 1 - 0.0 

- 1.1 

60.7 ± 5.3 

- 1 - 0.5 

- 3.5 

50.0 ± 7.3 

- 1 - 0.0 

- 13.0 

59.5 ± 11.6 

- 1 - 4.4 

- 0.0 


TABLE VII. Unpolarized cross sections in pb-GeV ^ with their statistical and asymmetric point-to-point 
systematic uncertainties for the Kin2 setting, for each bin in (j> (vertical) and —t (horizontal). 




4 > ( deg ) 

xb = 0.345 
= 2.218 GeV ^ 
t = - 0.176 GeV ^ 

Xb = 0.363 
= 2.318 GeV ^ 
t = - 0.232 GeV ^ 

Xb = 0.368 
= 2.348 GeV ^ 
t = - 0.279 GeV ^ 

Xb = 0.371 
= 2.360 GeV ^ 
t = - 0.325 GeV ^ 

Xb = 0.373 
= 2.375 GeV ^ 
t = - 0.372 GeV ^ 

7.5 

104.3 ± 3.6 

-t 1.1 

- 0.2 

80.3 ± 3.5 

- 1 - 0.0 

- 1.1 

66.6 ± 4.2 

-t 3.1 

- 0.0 

64.8 ± 5.8 

- 1 - 0.5 

- 1.9 

60.6 ± 8.4 

- 1 - 0.0 

- 8.2 

22.5 

108.6 ± 3.6 

-t 1.6 

- 0.9 

80.6 ± 3.3 

- 1 - 3.6 

- 0.0 

71.9 ± 4.0 

- 1 - 2.6 

- 0.6 

60.2 ± 4.9 

- 1 - 0.0 

- 3.9 

54.8 ± 6.3 

- 1 - 0.0 

- 3.3 

37.5 

100.3 ± 3.4 

- 1 - 2.6 

- 0.0 

70.5 ± 3.0 

-t 0.1 

- 1.8 

58.0 ± 3.2 

- 1 - 0.0 

- 5.8 

42.3 ± 3.4 

- 1 - 0.9 

- 0.6 

41.2 ± 4.1 

- 1 - 2.2 

- 1.4 

52.5 

91.2 ± 3.2 

-t 1.5 

- 1.0 

60.4 ± 2.7 

- 1 - 0.5 

- 1.3 

48.2 ± 2.8 

- 1 - 0.2 

- 0.7 

39.2 ± 3.0 

-t 3.1 

- 0.0 

38.8 ± 3.3 

-f 1.4 

- 0.0 

67.5 

76.8 ± 2.9 

- 1 - 0.0 

- 3.1 

51.3 ± 2.4 

- 1 - 0.6 

- 0.4 

42.2 ± 2.4 

- 1 - 0.2 

- 2.5 

36.8 ± 2.6 

- 1 - 3.0 

- 0.0 

29.5 ± 2.6 

-f 1.4 

- 0.1 

82.5 

73.0 ± 2.8 

- 1 - 0.4 

- 0.6 

46.3 ± 2.2 

-t 1.1 

- 0.1 

39.4 ± 2.2 

- 1 - 0.0 

- 3.0 

27.3 ± 2.1 

- 1 - 0.0 

- 1.0 

23.8 ± 2.2 

- 1 - 0.6 

- 0.2 

97.5 

56.5 ± 2.5 

- 1 - 0.4 

- 0.7 

42.6 ± 2.0 

- 1 - 0.5 

- 0.7 

29.5 ± 1.9 

-f 2.1 

- 0.2 

25.1 ± 1.9 

- 1 - 0.0 

- 0.8 

20.3 ± 1.9 

- 1 - 0.0 

- 1.1 

112.5 

55.7 ± 2.3 

- 1 - 0.6 

- 1.2 

35.3 ± 1.8 

- 1 - 0.8 

- 0.3 

28.7 ± 1.7 

- 1 - 0.0 

- 0.8 

24.3 ± 1.7 

- 1 - 0.4 

- 0.7 

16.6 ± 1.6 

- 1 - 0.0 

- 1.8 

127.5 

49.0 ± 2.2 

- 1 - 0.4 

- 1.1 

34.8 ± 1.7 

- 1 - 0.0 

- 1.5 

23.8 ± 1.5 

- 1 - 0.0 

- 0.9 

21.9 ± 1.6 

- 1 - 0.5 

- 0.4 

16.5 ± 1.5 

- 1 - 0.0 

- 2.6 

142.5 

47.0 ± 2.1 

-t 1.4 

- 0.0 

28.4 ± 1.6 

-t 0.1 

- 0.5 

24.8 ± 1.5 

- 1 - 0.0 

- 0.6 

20.3 ± 1.5 

- 1 - 0.0 

- 0.6 

14.9 ± 1.4 

- 1 - 0.0 

- 1.3 

157.5 

41.1 ± 2.0 

-t 0.1 

- 1.4 

30.9 ± 1.6 

-t 0.1 

- 0.9 

25.7 ± 1.6 

- 1 - 0.0 

- 1.4 

19.0 ± 1.5 

- 1 - 0.2 

- 0.2 

14.0 ± 1.5 

- 1 - 0.9 

- 0.0 

172.5 

41.6 ± 2.0 

- 1 - 0.0 

- 0.9 

25.7 ± 1.5 

-t 0.1 

- 0.6 

25.7 ± 1.6 

-t 0.1 

- 0.5 

20.4 ± 1.6 

- 1 - 0.9 

- 0.0 

13.7 ± 1.7 

- 1 - 0.0 

- 0.7 

187.5 

38.5 ± 1.9 

- 1 - 0.0 

- 1.1 

27.7 ± 1.5 

- 1 - 0.0 

- 0.9 

20.6 ± 1.5 

- 1 - 0.0 

- 0.6 

18.1 ± 1.5 

- 1 - 0.3 

- 0.6 

15.1 ± 1.7 

-f 1.1 

- 0.2 

202.5 

41.1 ± 1.9 

- 1 - 0.4 

- 1.1 

29.7 ± 1.6 

-t 0.1 

- 0.6 

20.7 ± 1.5 

-t 0.1 

- 0.4 

22.8 ± 1.6 

- 1 - 0.0 

- 1.7 

16.7 ± 1.6 

- 1 - 0.5 

- 0.3 

217.5 

45.3 ± 2.1 

- 1 - 0.3 

- 0.6 

31.4 ± 1.6 

-t 0.1 

- 0.6 

24.2 ± 1.5 

- 1 - 0.0 

- 1.5 

19.1 ± 1.5 

- 1 - 0.0 

- 1.0 

14.2 ± 1.4 

- 1 - 0.3 

- 0.5 

232.5 

51.3 ± 2.2 

- 1 - 0.0 

- 0.9 

34.4 ± 1.7 

- 1 - 0.3 

- 0.6 

28.3 ± 1.6 

- 1 - 0.0 

- 1.0 

21.9 ± 1.5 

- 1 - 0.0 

- 1.5 

19.0 ± 1.5 

- 1 - 0.0 

- 0.7 

247.5 

53.3 ± 2.3 

- 1 - 0.0 

- 1.6 

35.8 ± 1.8 

- 1 - 0.5 

- 0.3 

27.2 ± 1.7 

-t 1.3 

- 0.1 

20.8 ± 1.6 

- 1 - 0.0 

- 1.5 

18.5 ± 1.7 

-t 0.1 

- 2.4 

262.5 

62.2 ± 2.5 

- 1 - 0.0 

- 3.2 

42.4 ± 2.0 

-t 0.1 

- 1.1 

31.7 ± 1.8 

-t 0.1 

- 1.0 

25.6 ± 1.9 

- 1 - 0.3 

- 0.1 

19.6 ± 1.8 

- 1 - 0.0 

- 1.3 

277.5 

70.0 ± 2.7 

- 1 - 0.4 

- 1.0 

43.9 ± 2.1 

- 1 - 0.0 

- 1.6 

35.1 ± 2.0 

-t 0.1 

- 0.3 

31.5 ± 2.1 

- 1 - 0.0 

- 1.2 

24.2 ± 2.0 

- 1 - 0.0 

- 0.9 

292.5 

80.1 ± 2.9 

- 1 - 2.3 

- 0.0 

54.6 ± 2.3 

-t 1.0 

- 0.6 

41.8 ± 2.2 

-f 1.2 

- 0.1 

36.4 ± 2.4 

- 1 - 0.0 

- 1.9 

31.2 ± 2.5 

- 1 - 0.0 

- 1.4 

307.5 

87.4 ± 3.1 

-t 4.3 

- 0.0 

63.3 ± 2.6 

- 1 - 0.0 

- 3.3 

51.4 ± 2.7 

- 1 - 0.0 

- 2.6 

40.0 ± 2.7 

- 1 - 0.0 

- 1.8 

33.3 ± 2.9 

-t 1.6 

- 0.2 

322.5 

97.4 ± 3.3 

- 1 - 0.0 

- 1.8 

66.8 ± 2.9 

- 1 - 0.0 

- 3.3 

51.4 ± 2.9 

- 1 - 0.0 

- 4.9 

52.6 ± 3.4 

- 1 - 0.0 

- 2.1 

45.5 ± 4.0 

- 1 - 0.0 

- 7.6 

337.5 

99.2 ± 3.5 

- 1 - 2.7 

- 0.0 

75.2 ± 3.3 

- 1 - 0.0 

- 2.2 

56.8 ± 3.5 

- 1 - 0.0 

- 3.2 

55.7 ± 4.6 

- 1 - 0.0 

- 8.2 

55.2 ± 6.4 

- 1 - 0.7 

- 4.0 

352.5 

98.6 ± 3.6 

-t 1.7 

- 0.0 

77.3 ± 3.5 

- 1 - 0.0 

- 4.7 

60.0 ± 4.1 

- 1 - 2.0 

- 1.1 

59.1 ± 5.7 

- 1 - 3.2 

- 0.0 

55.4 ± 8.7 

-t 3.1 

- 1.3 


TABLE VIII. Unpolarized cross sections in pb-GeV with their statistical and asymmetric point-to-point 
systematic uncertainties for the Kin3 setting, for each bin in (j> (vertical) and —t (horizontal). 




if ) ( deg ) 

xb = 0.345 

q 2 = 1 453 GeV ^ 

t = - 0.170 GeV ^ 

Xb = 0.374 
= 1.552 GeV ^ 
t = - 0.232 GeV ^ 

Xb = 0.385 
= 1.589 GeV ^ 
t = - 0.278 GeV ^ 

Xb = 0.391 
= 1.608 GeV ^ 
t = - 0.323 GeV ^ 

Xb = 0.399 
= 1.633 GeV ^ 
t = - 0.370 GeV ^ 

7.5 

0.8 ± 7.9 

+ 1.0 

- 3.0 

5.7 ± 7.9 

- 1 - 0.0 

- 3.8 

- 7.2 ± 10.8 

- 1 - 0.0 

- 13.5 

4.2 ± 13.8 

- 1 - 4.5 

- 3.2 

12.7 ± 20.7 

+ 0.0 

- 16.2 

22.5 

5.2 ± 7.1 

+ 0.0 

- 9.4 

- 2.2 ± 7.1 

- 1 - 8.4 

- 0.0 

3.1 ± 8.9 

- 1 - 9.9 

- 1.8 

10.8 ± 13.1 

- 1 - 0.0 

- 11.6 

25.7 ± 16.3 

+ 2.6 

- 1.3 

37.5 

16.5 ± 6.3 

+ 0.0 

- 10.7 

4.3 ± 6.1 

-t 1.9 

- 3.4 

20.3 ± 6.9 

- 1 - 0.0 

- 5.1 

18.0 ± 8.7 

- 1 - 3.8 

- 0.0 

25.2 ± 11.2 

+ 0.0 

- 12.1 

52.5 

18.8 ± 5.2 

+ 1.1 

- 1.1 

13.4 ± 5.1 

- 1 - 3.2 

- 1.2 

18.3 ± 5.7 

- 1 - 0.0 

- 4.2 

18.0 ± 6.6 

- 1 - 7.8 

- 0.0 

16.3 ± 8.3 

+ 2.6 

- 1.8 

67.5 

19.9 ± 4.5 

+ 2.6 

- 0.0 

12.4 ± 4.6 

- 1 - 0.0 

- 4.0 

5.4 ± 4.8 

- 1 - 2.4 

- 0.0 

13.9 ± 5.4 

- 1 - 0.4 

- 4.6 

11.7 ± 6.2 

+ 0.0 

- 6.8 

82.5 

14.0 ± 4.0 

+ 1.6 

- 0.0 

12.2 ± 3.8 

-t 1.0 

- 0.0 

9.0 ± 4.1 

-t 6.1 

- 0.0 

5.7 ± 4.6 

- 1 - 3.6 

- 0.0 

6.9 ± 4.9 

+ 0.6 

- 2.2 

97.5 

15.8 ± 3.8 

+ 0.0 

- 2.8 

8.3 ± 3.5 

-t 1.6 

- 1.1 

8.6 ± 3.9 

-t 1.8 

- 2.9 

0.6 ± 4.2 

- 1 - 4.5 

- 0.0 

18.9 ± 4.7 

+ 1.4 

- 2.3 

112.5 

15.9 ± 3.6 

+ 1.9 

- 0.2 

13.6 ± 3.3 

- 1 - 0.7 

- 1.7 

- 2.2 ± 3.5 

-t 5.1 

- 0.0 

6.3 ± 3.9 

-t 1.3 

- 1.4 

2.2 ± 4.7 

+ 4.8 

- 0.8 

127.5 

12.7 ± 3.3 

+ 1.5 

- 0.3 

6.7 ± 3.1 

- 1 - 0.0 

- 4.9 

4.0 ± 3.2 

- 1 - 4.7 

- 0.0 

2.0 ± 3.9 

-t 0.1 

- 2.1 

7.1 ± 5.8 

+ 0.3 

- 7.2 

142.5 

10.0 ± 3.1 

+ 0.0 

- 3.0 

7.7 ± 2.8 

- 1 - 0.0 

- 3.4 

4.2 ± 3.3 

-t 1.8 

- 0.0 

1.2 ± 4.5 

- 1 - 0.6 

- 3.0 

- 1.5 ± 7.1 

+ 4.4 

- 2.8 

157.5 

6.8 ± 2.9 

+ 3.0 

- 0.0 

3.0 ± 2.7 

-f 1.1 

- 0.0 

8.0 ± 3.7 

- 1 - 0.9 

- 1.8 

5.9 ± 5.6 

- 1 - 4.0 

- 0.7 

9.1 ± 10.9 

+ 3.5 

- 4.0 

172.5 

- 3.5 ± 2.8 

+ 0.4 

- 1.6 

- 1.6 ± 2.7 

-t 1.5 

- 0.3 

- 2.0 ± 3.6 

- 1 - 5.0 

- 0.0 

- 1.8 ± 6.6 

- 1 - 2.3 

- 0.0 

- 18.3 ± 23.5 

+ 0.0 

- 34.1 

187.5 

- 1.9 ± 2.8 

+ 1.1 

- 0.3 

1.0 ± 2.6 

-t 1.0 

- 0.4 

0.2 ± 3.6 

-f 1.2 

- 0.8 

- 13.9 ± 6.1 

- 1 - 0.0 

- 4.2 

- 28.1 ± 15.5 

+ 5.6 

- 7.0 

202.5 

- 5.3 ± 2.9 

+ 1.1 

- 0.3 

- 1.8 ± 2.8 

- 1 - 4.0 

- 0.0 

0.9 ± 3.4 

- 1 - 0.0 

- 3.4 

1.8 ± 5.5 

- 1 - 2.0 

- 2.1 

0.6 ± 9.0 

+ 0.0 

- 8.1 

217.5 

- 4.4 ± 3.0 

+ 2.0 

- 0.8 

- 4.9 ± 2.9 

-f 1.2 

- 0.9 

- 1.5 ± 3.4 

-f 1.1 

- 1.9 

- 9.8 ± 4.8 

-f 1.2 

- 3.2 

0.1 ± 7.4 

+ 2.3 

- 0.8 

232.5 

- 8.0 ± 3.2 

+ 1.1 

- 0.4 

- 8.7 ± 2.9 

-f 1.1 

- 1.2 

- 5.2 ± 3.2 

-t 1.9 

- 0.0 

1.8 ± 3.9 

-t 1.8 

- 0.0 

- 1.9 ± 5.6 

+ 2.3 

- 0.2 

247.5 

- 11.8 ± 3.4 

+ 0.0 

- 5.7 

- 10.4 ± 3.2 

- 1 - 2.4 

- 0.0 

- 11.5 ± 3.5 

- 1 - 0.7 

- 0.2 

- 11.7 ± 4.1 

- 1 - 2.4 

- 0.0 

- 10.9 ± 5.1 

+ 0.8 

- 1.3 

262.5 

- 11.3 ± 3.6 

+ 0.0 

- 3.4 

- 7.9 ± 3.5 

- 1 - 2.4 

- 1.0 

- 6.2 ± 3.7 

- 1 - 0.8 

- 1.5 

- 15.0 ± 4.3 

- 1 - 0.4 

- 1.5 

- 5.5 ± 4.6 

+ 0.3 

- 5.1 

277.5 

- 24.0 ± 3.9 

+ 2.4 

- 0.2 

- 11.0 ± 3.6 

- 1 - 0.0 

- 3.3 

- 9.3 ± 3.9 

-t 1.0 

- 0.0 

- 12.1 ± 4.6 

-t 1.6 

- 2.1 

- 5.8 ± 4.9 

+ 0.9 

- 0.5 

292.5 

- 26.5 ± 4.4 

+ 2.4 

- 0.0 

- 16.4 ± 4.2 

- 1 - 0.9 

- 1.5 

- 7.5 ± 4.5 

- 1 - 0.9 

- 1.7 

- 4.2 ± 5.2 

- 1 - 0.0 

- 7.3 

- 11.7 ± 6.1 

+ 9.2 

- 0.0 

307.5 

- 16.6 ± 5.0 

+ 5.2 

- 1.5 

- 10.0 ± 4.7 

-t 3.1 

- 0.0 

- 2.5 ± 5.3 

- 1 - 0.0 

- 4.9 

1.8 ± 6.7 

- 1 - 0.0 

- 4.2 

- 11.7 ± 7.8 

+ 2.6 

- 0.9 

322.5 

- 18.7 ± 6.0 

+ 4.7 

- 0.0 

- 11.8 ± 6.0 

- 1 - 0.0 

- 5.5 

- 17.4 ± 6.9 

-t 7.5 

- 0.0 

- 15.0 ± 9.4 

- 1 - 7.9 

- 0.0 

- 14.9 ± 11.7 

+ 6.7 

- 0.5 

337.5 

- 5.7 ± 6.9 

+ 2.7 

- 2.3 

- 6.8 ± 7.0 

-t 4.3 

- 0.0 

- 3.0 ± 9.4 

- 1 - 0.0 

- 12.4 

- 8.7 ± 12.8 

- 1 - 9.5 

- 0.0 

- 28.7 ± 20.1 

- 1 - 16.1 

- 0.6 

352.5 

- 1.4 ± 7.6 

+ 1.6 

- 1.9 

- 4.3 ± 7.8 

+ 15.1 

- 0.0 

0.5 ± 11.2 

- 1 - 2.6 

- 3.8 

- 10.2 ± 15.5 

- 1 - 0.0 

- 7.2 

- 43.5 ± 23.1 

+ 8.2 

- 0.0 


TABLE IX. Cross-section differences for opposite beam helicities in pb-GeV”"^ with their statistical and 
asymmetric point-to-point systematic uncertainties for the Kinl setting, for each bin in (vertical) and —t 
(horizontal). 



<p ( deg ) 

xb = 0.343 
= 1.820 GeV ^ 
t = - 0.172 GeV ^ 

Xb = 0.368 

= 1.933 GeV ^ 

t = - 0.232 GeV ^ 

Xb = 0.375 

Q 2 = 1,964 GeV ^ 
t = - 0.278 GeV ^ 

Xb = 0.379 
= 1,986 GeV ^ 
t = - 0.323 GeV ^ 

Xb = 0.381 

= 1,999 GeV ^ 

t = - 0.371 GeV ^ 

7.5 

- 1.0 ± 6.0 

- 1.1 

3.1 ± 6.1 

- 4.8 

1.7 ± 7.7 + 

- 2.5 

7.3 ± 11.1 

- 0.0 

+ 14 

- 12.1 ± 16.0 

- 5.2 

22.5 

12.2 ± 5.7 

- 0.0 

- 1.9 ± 5.4 + 

- 1.4 

4.1 ± 6.2 

- 7.5 

- 0.1 ± 9.0 

- 3.7 

1.1 ± 12.0 

- 1.4 

37.5 

16.8 ± 4.9 

- 0.5 

12.0 ± 4.5 

- 2.0 

+ 42 

16.1 ± 5.2 ^ 

- 0.0 

7.9 ± 6.1 

- 3.4 

4.1 ± 7.5 

- 0.0 

52.5 

4 - 1 9 

15.3 ± 4.5 ^ 

- 0.9 

19.3 ± 4.2 

- 0.9 

15.0 ± 4.3 

- 0.0 

22.4 ± 5.1 

- 1.9 

+ 44 

7.3 ± 5.8 ^ 

- 0.2 

67.5 

15.8 ± 4.3 

- 1.2 

14.4 ± 3.7 

- 3.0 

1.8 ± 3.6 ^ 

- 1.1 

9.9 ± 4.0 

- 1.6 

9.9 ± 4.9 

- 1.1 

82.5 

20.5 ± 3.9 

- 0.0 

10.8 ± 3.2 

- 3.1 

+ 20 

10.9 ± 3.2 

- 0.0 

5.1 ± 3.5 

- 0.0 

10.8 ± 3.8 

- 1.4 

97.5 

15.1 ± 3.6 + 

- 0.0 

8.6 ± 2.8 

- 0.5 

4.8 ± 2.9 ^ 

- 1.7 

11.3 ± 3.0 

- 0.4 

+ 14 

4.7 ± 3.5 ^ 

- 0.3 

112.5 

6.7 ± 3.4 

- 1.7 

4.3 ± 2.7 + 

- 0.7 

8.6 ± 2.7 

- 1.3 

1.5 ± 2.7 

- 0.1 

+ 15 

10.5 ± 3.2 ^ 

- 1.2 

127.5 

14.9 ± 3.2 

- 1.1 

7.2 ± 2.4 + 

- 1.0 

6.4 ± 2.3 ^ 

- 0.0 

12.2 ± 2.6 

- 1.9 

2.1 ± 2.9 ■*“ 

- 0.0 

142.5 

6.9 ± 3.0 

- 0.0 

7.0 ± 2.3 

- 0.5 

6.9 ± 2.2 ^ 

- 0.0 

1.0 ± 2.5 

- 0.0 

3.7 ± 2.8 

- 3.6 

157.5 

-L 1 C; 

5.7 ± 3.0 ^ 

- 0.9 

- 1 - 1 3 

6.2 ± 2.3 

- 0.7 

5.4 ± 2.2 

- 1.1 

+ 14 

2.7 ± 2.6 

- 0.0 

0.1 ± 3.1 ■*“ 

- 1.6 

172.5 

0.7 ± 2.9 

- 1.0 

0.4 ± 2.2 

- 1.5 

- 1.5 ± 2.2 + 

- 1.4 

+ 21 

- 0.3 ± 2.5 ^ 

- 0.4 

- 7.0 ± 3.8 

- 0.0 

187.5 

±14 

- 0.3 ± 2.9 ^ 

- 0.0 

±14 

- 0.7 ± 2.2 ^ 

- 0.0 

- 1.4 ± 2.1 + 

- 0.5 

- 3.6 ± 2.7 

- 0.0 

2.5 ± 3.7 

- 0.0 

202.5 

- 6.5 ± 3.0 

- 0.7 

- 1.4 ± 2.3 

- 2.1 

- 0.5 ± 2.1 + 

- 1.3 

+ 14 

1.5 ± 2.6 ^ 

- 0.0 

0.6 ± 3.4 

- 1.3 

217.5 

- 8.1 ± 3.1 “*■ 

- 1.9 

- 2.6 ± 2.3 

- 2.2 

- 4.6 ± 2.2 ^ 

- 0.0 

2.3 ± 2.6 

- 0.8 

0.4 ± 3.1 

- 0.0 

232.5 

- 9.2 ± 3.1 

- 1.5 

- 6.7 ± 2.4 

- 2.1 

- 5.0 ± 2.3 ^ 

- 2.4 

0.0 ± 2.7 

- 0.0 

- 4.4 ± 2.7 

- 0.9 

247.5 

- 11.8 ± 3.2 

- 0.0 

- 10.6 ± 2.6 

- 1.1 

+ 04 

- 7.0 ± 2.8 ^ 

- 0.8 

- 8.1 ± 3.0 

- 1.6 

0.4 ± 3.1 

- 2.5 

262.5 

- 17.5 ± 3.4 + 3-° 
- 0.0 

4 - 1 9 

- 9.2 ± 2.8 

- 0.0 

- 11.0 ± 3.0 ^ 

- 1.0 

- 7.1 ± 3.1 

- 0.0 

- 2.5 ± 3.2 

- 0.7 

277.5 

- 16.0 ± 3.8 

- 1.4 

- 15.5 ± 3.0 

- 3.5 

- 11.9 ± 3.2 + 

- 0.0 

0.6 ± 3.5 

- 0.0 

- 6.1 ± 3.6 

- 2.1 

292.5 

- 17.2 ± 4.0 

- 0.4 

- 16.0 ± 3.5 

- 0.0 

- 9.9 ± 3.5 + 

- 1.3 

- 8.9 ± 4.0 

- 1.5 

- 1.6 ± 4.4 

- 1.1 

307.5 

±17 

- 16.3 ± 4.5 ^ 

- 0.6 

- 10.3 ± 3.9 

- 0.0 

- 6.0 ± 4.1 +°-° 
- 3.9 

- 8.4 ± 5.1 

- 0.0 

- 7.1 ± 5.4 

- 1.9 

322.5 

±07 

- 11.6 ± 4.9 

- 1.3 

- 16.5 ± 4.5 

- 0.0 

- 9.6 ± 4.8 ^ 

- 0.0 

- 20.1 ± 6.3 

- 3.6 

- 7.9 ± 7.6 

- 3.5 

337.5 

- 1.9 ± 5.5 

- 1.0 

- 7.2 ± 5.5 

- 2.3 

- 17.2 ± 6.2 ^ 

- 8.0 

- 3.3 ± 8.8 

- 6.7 

±17 

- 1.2 ± 10.9 

- 15.4 

352.5 

- 3.6 ± 5.9 + 

- 2.2 

- 0.7 ± 6.1 

- 1.1 

1.8 ± 7.4 

- 1.3 

+ 30 

- 2.6 ± 10.6 

- 8.5 

- 5.1 ± 16.7 + 

- 0.0 


TABLE X. Cross-section differences for opposite beam helicities in pb-GeV”"* with their statistical and 
asymmetric point-to-point systematic uncertainties for the Kin2 setting, for each bin in (vertical) and —t 
(horizontal). 




if ) ( deg ) 

xb = 0.345 
= 2.218 GeV ^ 
t = - 0.176 GeV ^ 

Xb = 0.363 
= 2.318 GeV ^ 
t = - 0.232 GeV ^ 

Xb = 0.368 

Q 2 = 2.348 GeV ^ 
t = - 0.279 GeV ^ 

Xb = 0.371 
= 2.360 GeV ^ 
t = - 0.325 GeV ^ 

Xb = 0.373 
= 2.375 GeV ^ 
t = - 0.372 GeV ^ 

7.5 

5.2 ± 5.1 

- 5.4 

7.4 ± 4.9 

- 1.8 

11.5 ± 5.9 

- 0.1 

1.7 ± 7.3 

- 0.0 

9.4 ± 10.9 

- 2.5 

22.5 

1.8 ± 5.1 

- 0.0 

15.4 ± 4.6 + 

- 1.3 

8.6 ± 5.3 + 

- 0.4 

4.4 ± 6.2 

- 4.3 

6.4 ± 8.2 

- 1.9 

37.5 

8.0 ± 4.7 + 

- 0.7 

8.3 ± 4.1 

- 0.0 

15.4 ± 4.2 + 

- 1.9 

1.9 ± 4.4 

- 5.3 

10.0 ± 5.3 

- 3.6 

52.5 

^57 

16.7 ± 4.5 ^ 

- 0.1 

±15 

11.5 ± 3.6 ^ 

- 1.5 

± 1 5 

10.5 ± 3.6 ^ 

- 0.2 

9.6 ± 3.8 

- 2.3 

21.8 ± 4.3 

- 0.4 

67.5 

21.1 ± 4.1 

- 0.0 

±37 

17.1 ± 3.2 ^ 

- 0.0 

± 1 2 

16.7 ± 3.2 ^ 

- 0.0 

15.0 ± 3.4 + 

- 1.4 

6.2 ± 3.3 

- 3.5 

82.5 

16.0 ± 4.0 ^ 

- 0.1 

6.6 ± 2.9 

- 0.0 

11.9 ± 2.8 

- 0.4 

±17 

9.3 ± 2.6 

- 0.2 

9.3 ± 2.7 + 

- 0.5 

97.5 

+ 1 9 

10.1 ± 3.5 ^ 

- 0.1 

12.3 ± 2.7 ■*“ 

- 0.7 

±14 

9.8 ± 2.4 ^ 

- 0.0 

9.9 ± 2.4 

- 0.0 

9.0 ± 2.4 

- 1.8 

112.5 

20.7 ± 3.4 ^ 

- 1.8 

11.6 ± 2.4 + 

- 1.3 

8.0 ± 2.2 

- 0.9 

3.7 ± 2.2 

- 0.3 

- 1.1 ± 2.0 + 

- 0.0 

127.5 

9.3 ± 3.2 + 

- 1.9 

±07 

6.4 ± 2.3 

- 1.4 

3.5 ± 1.9 

- 1.0 

±04 

6.6 ± 2.0 ^ 

- 0.0 

1.6 ± 1.8 

- 0.4 

142.5 

4.8 ± 3.1 + 

- 0.1 

5.5 ± 2.1 + 

- 0.2 

3.7 ± 2.0 

- 0.0 

3.6 ± 1.9 

- 1.5 

3.2 ± 1.8 

- 0.0 

157.5 

±14 

- 0.4 ± 3.0 ^ 

- 0.0 

2.7 ± 2.2 

- 0.1 

4.0 ± 2.0 

- 0.5 

5.5 ± 1.8 

- 0.2 

±07 

3.6 ± 1.9 ^ 

- 0.0 

172.5 

±29 

3.3 ± 3.0 ^ 

- 0.2 

0.0 ± 2.0 

- 0.0 

±07 

- 0.3 ± 2.0 

- 0.1 

2.7 ± 1.9 

- 0.9 

2.1 ± 2.0 

- 0.2 

187.5 

- 0.6 ± 2.8 ^ 

- 1.4 

0.2 ± 2.1 

- 0.2 

1.0 ± 1.9 

- 1.2 

- 0.6 ± 1.9 

- 0.5 

4.4 ± 2.0 

- 1.0 

202.5 

- 2.6 ± 2.9 ^ 

- 2.0 

±04 

- 0.7 ± 2.1 ^ 

- 0.7 

- 1.8 ± 1.8 

- 0.5 

- 4.2 ± 2.0 

- 0.0 

- 0.1 ± 2.0 

- 2.5 

217.5 

- 7.5 ± 3.0 + 

- 0.0 

- 5.4 ± 2.2 + 

- 0.5 

- 5.1 ± 2.0 + 

- 0.0 

- 5.2 ± 1.8 

- 0.7 

- 3.5 ± 1.8 

- 0.1 

232.5 

- 13.1 ± 3.2 + 

- 0.0 

- 4.2 ± 2.3 

- 1.2 

- 4.7 ± 2.1 

- 0.8 

±14 

- 5.2 ± 1.9 ^ 

- 0.0 

- 4.0 ± 1.9 

- 0.7 

247.5 

±17 

- 6.9 ± 3.3 

- 0.0 

- 9.5 ± 2.4 + 

- 0.7 

- 3.8 ± 2.1 

- 0.0 

0.6 ± 2.1 

- 0.4 

- 1.8 ± 2.1 

- 1.0 

262.5 

± 2 1 

- 17.5 ± 3.6 ^ 

- 0.0 

- 9.0 ± 2.7 

- 1.4 

- 10.5 ± 2.4 

- 3.0 

- 6.8 ± 2.4 

- 1.7 

- 3.4 ± 2.3 

- 1.1 

277.5 

±14 

- 13.0 ± 3.8 

- 0.7 

±05 

- 16.0 ± 2.8 

- 1.2 

- 9.8 ± 2.6 

- 0.0 

- 14.2 ± 2.7 + 

- 1.5 

±07 

- 5.8 ± 2.6 ^ 

- 0.3 

292.5 

±29 

- 13.2 ± 4.2 ^ 

- 0.6 

±04 

- 12.6 ± 3.1 

- 1.8 

- 11.0 ± 2.9 

- 0.9 

- 9.6 ± 3.1 

- 3.6 

- 11.3 ± 3.1 

- 1.8 

307.5 

±04 

- 16.6 ± 4.3 ^ 

- 2.1 

±44 

- 14.0 ± 3.5 ^ 

- 0.0 

±17 

- 15.8 ± 3.5 

- 0.5 

±07 

- 10.6 ± 3.6 

- 0.6 

- 8.8 ± 3.7 

- 0.5 

322.5 

- 15.2 ± 4.6 + 

- 0.0 

- 14.0 ± 3.8 

- 1.5 

- 9.4 ± 3.8 

- 0.0 

±44 

- 14.2 ± 4.5 

- 0.0 

- 5.2 ± 5.1 

- 2.6 

337.5 

- 3.1 ± 4.8 + 

- 0.4 

1,4 ± 4.4 + 

- 0.8 

- 5.0 ± 4.7 + 

- 1.7 

±07 

- 0.5 ± 5.9 ^ 

- 2.6 

- 8.1 ± 8.1 + 

- 4.3 

352.5 

- 1.2 ± 4.9 + 

- 1.5 

- 1.1 ± 4.8 

- 0.0 

2.7 ± 5.5 

- 3.8 

- 0.5 ± 7.3 + 

- 1.5 

7.3 ± 11.5 

- 1.1 


TABLE XI. Cross-section differences for opposite beam helicities in pb-GeV”"^ with their statistical and 
asymmetric point-to-point systematic uncertainties for the Kin3 setting, for each bin in (vertical) and —t 
(horizontal). 




0 ( deg ) 

xb = 0.378 
= 2.012 GeV ^ 
t = - 0.192 GeV ^ 

Xb = 0.392 
= 2.054 GeV ^ 
t = - 0.233 GeV ^ 

xg = 0.398 

Q2 = 2.074 GeV ^ 
t = - 0.279 GeV ^ 

Xg = 0.400 
= 2.084 GeV ^ 
t = - 0.324 GeV ^ 

Xg = 0.401 
= 2.091 GeV ^ 
t = - 0.371 GeV ^ 

7.5 

64.0 ± 6.3 

-t 1.1 

- 0.9 

57.0 ± 3.8 

-t 0.1 

- 1.2 

56.3 ± 4.4 

-f 1.4 

- 1.7 

50.5 ± 5.9 

-f 2.1 

- 1.3 

53.3 ± 9.0 

- 1 - 8.9 

- 0.0 

22.5 

68.2 ± 6.2 

-t 1.8 

- 2.0 

55.1 ± 3.6 

- 1 - 0.0 

- 2.5 

45.2 ± 3.8 

-t 0.1 

- 0.9 

51.7 ± 5.1 

- 1 - 0.8 

- 1.7 

52.0 ± 6.9 

- 1 - 3.4 

- 0.3 

37.5 

52.8 ± 5.6 

- 1 - 0.0 

- 4.6 

48.2 ± 3.3 

- 1 - 0.5 

- 1.3 

44.9 ± 3.5 

-f 1.4 

- 0.5 

36.2 ± 4.0 

-t 1.0 

- 0.0 

36.9 ± 4.8 

- 1 - 0.0 

- 3.8 

52.5 

52.4 ± 5.8 

- 1 - 0.6 

- 2.8 

51.4 ± 3.3 

- 1 - 3.3 

- 0.1 

41.2 ± 3.4 

- 1 - 0.0 

- 4.4 

36.4 ± 3.8 

- 1 - 0.0 

- 2.7 

33.8 ± 4.1 

- 1 - 0.0 

- 4.5 

67.5 

53.0 ± 5.4 

-t 1.3 

- 0.5 

42.0 ± 3.0 

- 1 - 0.0 

- 3.0 

34.6 ± 3.0 

- 1 - 0.0 

- 1.5 

30.7 ± 3.2 

- 1 - 0.0 

- 3.2 

31.0 ± 3.5 

- 1 - 0.0 

- 2.0 

82.5 

47.3 ± 5.1 

- 1 - 2.3 

- 0.0 

41.4 ± 2.9 

- 1 - 0.0 

- 3.0 

30.9 ± 2.9 

- 1 - 3.3 

- 0.0 

24.9 ± 3.0 

-t 1.8 

- 0.0 

22.0 ± 3.0 

- 1 - 0.0 

- 3.9 

97.5 

48.2 ± 4.9 

- 1 - 3.3 

- 0.0 

32.8 ± 2.7 

- 1 - 0.3 

- 2.6 

30.5 ± 2.7 

- 1 - 0.5 

- 2.0 

26.4 ± 2.9 

- 1 - 0.4 

- 1.4 

27.4 ± 3.2 

- 1 - 0.0 

- 2.7 

112.5 

46.9 ± 4.7 

- 1 - 4.8 

- 1.3 

33.9 ± 2.6 

- 1 - 0.0 

- 2.9 

28.9 ± 2.5 

- 1 - 0.0 

- 1.2 

20.0 ± 2.5 

- 1 - 0.0 

- 1.2 

24.0 ± 3.1 

-t 1.3 

- 0.2 

127.5 

45.8 ± 4.6 

- 1 - 0.0 

- 1.5 

28.0 ± 2.3 

- 1 - 0.0 

- 2.3 

23.1 ± 2.2 

- 1 - 0.0 

- 3.4 

24.5 ± 2.6 

- 1 - 0.0 

- 1.8 

19.7 ± 3.0 

- 1 - 0.0 

- 1.5 

142.5 

36.5 ± 4.2 

-t 1.0 

- 0.0 

27.3 ± 2.3 

- 1 - 0.0 

- 1.0 

23.2 ± 2.2 

- 1 - 0.0 

- 1.4 

18.9 ± 2.4 

- 1 - 0.3 

- 0.7 

16.2 ± 3.2 

- 1 - 0.0 

- 2.6 

157.5 

33.2 ± 4.0 

- 1 - 0.7 

- 1.3 

27.6 ± 2.2 

- 1 - 0.0 

- 1.7 

25.4 ± 2.3 

- 1 - 0.0 

- 2.8 

20.0 ± 2.7 

- 1 - 0.0 

- 1.6 

15.5 ± 3.6 

- 1 - 2.7 

- 0.0 

172.5 

36.3 ± 4.1 

- 1 - 0.0 

- 3.1 

26.7 ± 2.2 

-t 1.6 

- 0.3 

24.2 ± 2.2 

- 1 - 0.0 

- 2.3 

18.5 ± 2.6 

- 1 - 0.3 

- 1.2 

15.9 ± 3.9 

-t 0.1 

- 1.1 

187.5 

35.6 ± 4.1 

- 1 - 0.3 

- 0.8 

26.6 ± 2.2 

- 1 - 0.2 

- 1.0 

22.4 ± 2.1 

- 1 - 0.0 

- 0.9 

16.6 ± 2.6 

- 1 - 0.5 

- 0.4 

17.4 ± 4.0 

- 1 - 0.0 

- 3.5 

202.5 

37.7 ± 4.3 

- 1 - 0.0 

- 4.5 

30.9 ± 2.3 

-t 1.3 

- 0.0 

20.4 ± 2.0 

- 1 - 0.0 

- 1.1 

22.8 ± 2.8 

-t 1.5 

- 0.3 

17.2 ± 3.4 

- 1 - 0.0 

- 3.0 

217.5 

43.6 ± 4.5 

-t 1.1 

- 1.0 

28.3 ± 2.3 

- 1 - 0.8 

- 0.0 

29.1 ± 2.4 

- 1 - 0.0 

- 5.0 

22.1 ± 2.6 

- 1 - 0.0 

- 2.9 

22.2 ± 3.2 

- 1 - 0.3 

- 0.9 

232.5 

38.1 ± 4.4 

- 1 - 0.0 

- 3.5 

25.6 ± 2.3 

- 1 - 0.0 

- 1.5 

21.6 ± 2.2 

- 1 - 0.0 

- 1.2 

20.8 ± 2.6 

-t 1.0 

- 0.0 

17.1 ± 2.7 

-t 1.5 

- 0.1 

247.5 

37.7 ± 4.3 

- 1 - 0.9 

- 1.8 

30.5 ± 2.5 

-t 1.7 

- 0.0 

26.7 ± 2.5 

- 1 - 0.0 

- 3.1 

28.1 ± 2.9 

- 1 - 0.0 

- 4.2 

22.6 ± 3.0 

- 1 - 0.0 

- 4.7 

262.5 

43.4 ± 4.6 

- 1 - 0.0 

- 8.3 

34.3 ± 2.7 

- 1 - 0.7 

- 1.1 

30.0 ± 2.7 

- 1 - 0.0 

- 2.2 

24.8 ± 2.9 

-f 1.2 

- 0.2 

19.5 ± 2.9 

- 1 - 0.4 

- 1.6 

277.5 

44.1 ± 4.9 

-t 1.0 

- 0.8 

34.5 ± 2.8 

- 1 - 0.0 

- 0.9 

34.9 ± 2.8 

- 1 - 0.0 

- 4.2 

28.2 ± 3.1 

- 1 - 0.4 

- 1.1 

23.7 ± 3.1 

- 1 - 0.3 

- 1.0 

292.5 

50.5 ± 5.3 

- 1 - 2.3 

- 0.0 

41.8 ± 3.0 

- 1 - 0.2 

- 0.6 

30.5 ± 2.7 

-f 1.4 

- 0.2 

27.0 ± 3.0 

- 1 - 3.6 

- 0.0 

27.2 ± 3.5 

- 1 - 0.4 

- 1.4 

307.5 

57.7 ± 5.8 

- 1 - 0.3 

- 3.4 

43.2 ± 3.2 

- 1 - 0.0 

- 3.8 

35.9 ± 3.0 

-t 1.6 

- 0.0 

39.3 ± 3.6 

- 1 - 0.0 

- 1.4 

29.7 ± 3.8 

- 1 - 0.2 

- 2.8 

322.5 

62.1 ± 5.9 

- 1 - 0.9 

- 0.9 

53.2 ± 3.5 

- 1 - 0.0 

- 2.5 

41.0 ± 3.3 

- 1 - 0.0 

- 0.8 

39.5 ± 4.1 

- 1 - 0.6 

- 2.2 

31.6 ± 4.8 

- 1 - 0.0 

- 4.5 

337.5 

62.3 ± 5.9 

-t 1.5 

- 0.0 

53.2 ± 3.5 

- 1 - 2.3 

- 0.3 

46.9 ± 3.8 

-f 1.4 

- 0.7 

41.4 ± 5.0 

- 1 - 0.9 

- 1.2 

33.4 ± 6.4 

- 1 - 4.6 

- 0.0 

352.5 

65.4 ± 6.0 

- 1 - 0.0 

- 1.5 

55.5 ± 3.7 

- 1 - 0.0 

- 1.5 

44.0 ± 4.1 

- 1 - 0.3 

- 1.6 

37.4 ± 5.5 

- 1 - 0.0 

- 9.8 

43.6 ± 8.8 

- 1 - 2.9 

- 0.0 


TABLE XII. Unpolarized cross sections in pb-GeV ^ with their statistical and asymmetric point-to-point 
systematic uncertainties for the KinX2 setting, for each bin in (j) (vertical) and —t (horizontal). 




(/> ( deg ) 

xg = 0.336 
= 2.161 GeV ^ 
t = - 0.171 GeV ^ 

xq = 0.342 
= 2.190 GeV ^ 
t = - 0.231 GeV ^ 

xg = 0.343 

Q 2 = 2.194 GeV ^ 
t = - 0.278 GeV ^ 

xq = 0.342 
= 2.191 GeV ^ 
t = - 0.324 GeV ^ 

xq = 0.342 

Q 2 = 2.193 GeV ^ 
t = - 0.371 GeV ^ 

7.5 

120.7 ± 4.7 

- 0.4 

4 - 2 1 

100.1 ± 7.6 ^ 

- 2.5 

- 

- 

- 

22.5 

- 1-3 7 

132.5 ± 4.9 

- 0.0 

±45 

106.9 ± 6.8 

- 0.5 

- 

- 

- 

37.5 

118.7 ± 4.5 

- 0.0 

4-0 2 

93.9 ± 6.1 ^ 

- 1.6 

58.7 ± 7.0 

- 4.3 

4-0 7 

53.8 ± 10.1 ^ 

- 9.3 

4-0 9 

68.8 ± 21.3 

- 14.6 

52.5 

4 - 1 3 

99.8 ± 4.1 ^ 

- 1.9 

76.9 ± 5.0 

- 2.7 

59.3 ± 5.8 

- 1.1 

48.6 ± 7.3 

- 0.3 

37.1 ± 10.4 

- 9.7 

67.5 

88.4 ± 3.8 

- 5.5 

67.7 ± 4.5 

- 1.3 

P , ± b .9 

54.4 ± 5.1 

- 0.9 

55.8 ± 6.2 

- 0.0 

40.6 ± 7.0 

- 0.0 

82.5 

81.4 ± 3.6 

- 0.7 

55.1 ± 3.9 

- 0.0 

47.8 ± 4.2 

- 2.6 

±07 

34.2 ± 4.4 ^ 

- 1.8 

4 - 2 1 

29.0 ± 5.0 ■ 

- 1.0 

97.5 

-L 1 c ; 

62.7 ± 3.1 ^ 

- 0.8 

4-12 

53.8 ± 3.7 ^ 

- 0.7 

36.7 ± 3.6 

- 0.0 

30.9 ± 3.8 

- 1.1 

27.0 ± 4.3 

- 2.9 

112.5 

61.3 ± 3.0 

- 1.5 

±07 

40.3 ± 3.1 ^ 

- 0.7 

36.5 ± 3.3 

- 1.4 

28.6 ± 3.4 

- 1.2 

21.0 ± 3.5 

- 3.1 

127.5 

58.1 ± 3.0 + 

- 1.9 

37.3 ± 2.8 

- 3.7 

25.9 ± 2.6 

- 1.2 

27.6 ± 3.1 + 

- 0.6 

18.8 ± 2.9 

- 3.4 

142.5 

52.9 ± 2.7 

- 0.0 

33.6 ± 2.7 

- 1.7 

27.3 ± 2.7 

- 0.5 

26.7 ± 3.0 

- 2.0 

17.9 ± 2.8 

- 0.9 

157.5 

43.9 ± 2.6 

- 0.8 

±07 

34.4 ± 2.8 

- 0.6 

29.2 ± 2.8 

- 2.8 

24.8 ± 2.8 

- 0.0 

18.7 ± 3.0 

- 0.0 

172.5 

45.4 ± 2.5 + 

- 1.2 

4-0 2 

31.9 ± 2.6 ^ 

- 0.6 

4-0 3 

29.9 ± 2.8 

- 0.2 

4-17 

23.5 ± 2.8 

- 0.0 

17.2 ± 3.2 

- 0.3 

187.5 

41.0 ± 2.4 

- 0.1 

4-0 5 

32.6 ± 2.7 

- 0.3 

27.7 ± 2.8 

- 2.1 

±14 

21.6 ± 2.8 

- 0.0 

16.2 ± 2.9 

- 0.0 

202.5 

44.2 ± 2.5 

- 0.9 

37.1 ± 2.8 

- 1.6 

23.7 ± 2.6 

- 0.5 

29.2 ± 3.0 

- 3.1 

4 - 2 1 

21.3 ± 3.0 

- 0.5 

217.5 

51.6 ± 2.7 

- 0.8 

39.5 ± 2.9 

- 2.0 

31.2 ± 2.8 

- 0.7 

23.0 ± 2.7 

- 2.1 

17.7 ± 2.8 

- 0.8 

232.5 

57.7 ± 2.9 

- 0.8 

41.5 ± 3.0 

- 0.0 

33.4 ± 2.9 

- 1.9 

29.3 ± 3.0 

- 0.9 

24.2 ± 3.0 

- 1.3 

247.5 

56.9 ± 2.9 

- 2.3 

±14 

43.0 ± 3.1 ^ 

- 0.0 

30.9 ± 3.0 

- 0.4 

21.4 ± 3.1 

- 0.7 

±04 

22.5 ± 3.4 ^ 

- 3.9 

262.5 

69.3 ± 3.2 

- 3.7 

50.9 ± 3.5 

- 0.8 

4-0 7 

36.4 ± 3.4 ^ 

- 1.2 

30.4 ± 3.9 

- 0.3 

18.1 ± 3.6 

- 2.6 

277.5 

±07 

77.0 ± 3.5 ^ 

- 0.4 

4-0 5 

54.2 ± 3.7 

- 2.0 

±04 

48.5 ± 3.9 ^ 

- 1.8 

43.6 ± 4.4 + 

- 3.2 

32.8 ± 4.5 

- 2.5 

292.5 

95.5 ± 3.9 + 

- 0.0 

67.8 ± 4.2 

- 1.4 

51.6 ± 4.3 

- 0.2 

55.0 ± 5.6 

- 1.2 

35.5 ± 6.3 

- 1.9 

307.5 

±48 

103.1 ± 4.0 ^ 

- 0.0 

84.2 ± 4.8 ■*“ 

- 6.6 

60.3 ± 5.4 

- 3.9 

39.5 ± 6.2 

- 2.9 

±41 

43.1 ± 8.7 ^ 

- 0.3 

322.5 

±04 

114.1 ± 4.4 ^ 

- 1.7 

90.7 ± 5.5 

- 6.5 

62.2 ± 6.4 

- 8.6 

70.1 ± 9.9 

- 7.6 

49.6 ± 15.5 

- 10.8 

337.5 

117.8 ± 4.6 

- 0.6 

83.8 ± 6.6 

- 5.6 

- 

- 

- 

352.5 

4-0 2 

120.3 ± 4.8 

- 0.7 

87.9 ± 7.4 

- 10.9 

- 

- 

- 


TABLE XIII. Unpolarized cross sections in pb-GeV with their statistical and asymmetric point-to-point 
systematic uncertainties for the KinX3 setting, for each bin in (j) (vertical) and —t (horizontal). 



0 ( deg ) 

xb = 0.378 
= 2.012 GeV ^ 
t = - 0.192 GeV ^ 

xg = 0.392 
= 2.054 GeV ^ 
t = - 0.233 GeV ^ 

Xg = 0.398 

Q 2 = 2.074 GeV ^ 
t = - 0.279 GeV ^ 

Xg = 0.400 
= 2.084 GeV ^ 
t = - 0.324 GeV ^ 

Xg = 0.401 

Q 2 = 2.091 GeV ^ 
t = - 0.371 GeV ^ 

7.5 

- 2.8 ± 9.1 + 

- 2.5 

- 2.4 ± 5.1 + 

- 3.5 

- 0.2 ± 5.9 + 

- 1.0 

3.8 ± 7.7 + 

- 0.0 

- 9.0 ± 11.9 + 

- 3.2 

22.5 

+ 07 

2.7 ± 8.8 ^ 

- 5.6 

- 0.6 ± 4.8 + 

- 2.4 

1.4 ± 5.0 +°-° 
- 4.0 

- 1.2 ± 6.8 + 

- 3.5 

+ 27 

0.9 ± 9.2 ^ 

- 1.1 

37.5 

- 6.0 ± 7.9 + 

- 0.0 

4.1 ± 4.4 + 

- 1.1 

13.3 ± 4.6 + 

- 0.0 

5.5 ± 5.2 + 2-0 
- 3.0 

+ 14 

0.4 ± 6.2 ^ 

- 0.0 

52.5 

7.2 ± 8.1 

- 2.5 

13.7 ± 4.5 + 2-6 
- 0.4 

11.2 ± 4.4 + 2-3 
- 0.0 

+ 19 

12.3 ± 4.8 

- 1.2 

8.3 ± 5.3 + 2-0 
- 0.4 

67.5 

- 5.9 ± 7.7 + 2.2 
- 0.7 

5.0 ± 4.1 + ^ + 
- 0.0 

6.0 ± 3.9 + 

- 1.4 

8.4 ± 4.2 +°-'^ 
- 2.0 

5.9 ± 4.5 +°-® 
- 1.0 

82.5 

+ 07 

18.2 ± 7.2 ^ 

- 1.9 

12.4 ± 4.0 + 

- 3.3 

10.6 ± 3.7 + °-2 
- 1.8 

3.5 ± 3.7 + 

- 1.4 

7.9 ± 3.8 +°-° 
- 1.7 

97.5 

- 0.8 ± 7.3 + 

- 0.6 

4.6 ± 3.6 + ^'2 
- 0.0 

4.2 ± 3.5 + °-2 
- 1.8 

7.2 ± 3.7 + 

- 0.0 

3.3 ± 4.0 + 

- 0.2 

112.5 

14.0 ± 6.7 + 2-8 
- 0.0 

- 0.6 ± 3.6 + 

- 2.5 

7.2 ± 3.3 + 

- 1.1 

- 1.1 ± 3.1 + 2 -^ 
- 0.0 

6.5 ± 3.9 +^-® 
- 0.5 

127.5 

- 0.3 ± 6.7 + 2 -^ 
- 0.7 

3.0 ± 3.1 + 

- 1.6 

7.3 ± 2.8 + 

- 0.0 

11.1 ± 3.2 +°-° 
- 4.8 

- 1.5 ± 3.6 + 2-9 
- 0.0 

142.5 

1.6 ± 6.0 + 

- 4.0 

4.8 ± 3.1 + °-2 
- 1.7 

4.8 ± 2.8 + 2 -® 
- 0.0 

3.3 ± 3.0 + 

- 0.0 

6.1 ± 3.8 + 

- 3.8 

157.5 

- 9.0 ± 5.9 + 2-5 
- 2.4 

4.9 ± 3.0 + 

- 1.2 

5.6 ± 2.9 + 2 + 
- 0.8 

5.3 ± 3.3 +°-^ 
- 1.5 

- 2.7 ± 4.2 +°-® 
- 0.6 

172.5 

- 4.9 ± 5.9 + 

- 3.1 

1.0 ± 3.0 + 

- 1.5 

0.6 ± 2.8 + +3 
- 0.1 

4 - 1 ^ 

- 0.7 ± 3.2 ^ ° 

- 0.3 

- 8.4 ± 4.5 +°-^ 
- 2.0 

187.5 

+ 14 

- 2.7 ± 6.0 ^ 

- 3.3 

2.5 ± 3.0 + 

- 1.1 

1.4 ± 2.7 +°-^ 
- 1.1 

- 0.5 ± 3.2 + 2-2 
- 0.7 

+ 47 

2.9 ± 4.7 ^ 

- 0.0 

202.5 

6.7 ± 6.3 + 

- 4.1 

- 1.0 ± 3.2 + 

- 3.2 

1.7 ± 2.6 + 

- 1.4 

1.0 ± 3.4 + +® 
- 0.0 

3.8 ± 4.1 + 

- 1.4 

217.5 

- 0.8 ± 6.6 + 

- 10.0 

1.3 ± 3.1 + 

- 2.6 

- 5.0 ± 3.1 + 

- 0.0 

- 1.2 ± 3.2 + 

- 1.6 

0.8 ± 3.9 + 2 -® 
- 0.0 

232.5 

- 8.7 ± 6.3 + 2 -® 
- 0.0 

- 1.1 ± 3.1 + 

- 4.4 

- 4.3 ± 2.8 + 

- 1.0 

0.3 ± 3.2 + ^'2 
- 0.7 

- 6.6 ± 3.2 + 

- 0.2 

247.5 

- 12.9 ± 6.2 + 

- 4.0 

- 10.3 ± 3.4 + 

- 2.8 

+ 04 

- 6.8 ± 3.2 ^ 

- 0.2 

- 6.6 ± 3.7 + 

- 1.2 

- 3.0 ± 3.7 +°-° 
- 3.2 

262.5 

- 4.0 ± 6.5 + 

- 3.8 

- 2.4 ± 3.7 + 2 -® 
- 0.0 

+ 1 9 

- 8.4 ± 3.5 ^ 

- 0.4 

+ 1 3 

- 7.9 ± 3.6 ^ 

- 1.1 

- 0.1 ± 3.7 + 

- 0.4 

277.5 

- 1.2 ± 6.9 + 

- 2.9 

- 12.5 ± 3.7 + 

- 3.0 

- 9.4 ± 3.7 + 2 -^ 
- 0.6 

4.4 ± 3.8 +^-'^ 
- 0.0 

+ 04 

- 3.2 ± 3.9 ^ 

- 1.0 

292.5 

- 6.0 ± 7.4 + 

- 1.6 

- 13.8 ± 4.1 + 2 -® 
- 0.0 

- 7.3 ± 3.5 + +3 
- 0.3 

- 4.5 ± 3.9 + 

- 0.7 

- 7,4 ± 4.4 + 2 -"^ 
- 0.6 

307.5 

+ 17 

0.6 ± 8.3 

- 5.6 

- 5.3 ± 4.2 + 

- 0.0 

- 6.2 ± 4.0 + 

- 1.8 

+ 47 

- 4.1 ± 4.8 ^ 

- 0.2 

- 7.2 ± 4.9 +°-® 
- 1.7 

322.5 

- 2.6 ± 8.3 + 

- 8.0 

- 15.7 ± 4.6 + 

- 1.2 

- 10.4 ± 4.3 + 

- 1.1 

- 14.6 ± 5.3 + 

- 3.9 

- 5.8 ± 6.2 +^-° 
- 2.7 

337.5 

- 14.3 ± 8.6 + 

- 1.6 

- 1.8 ± 4.8 + 

- 0.8 

- 14.9 ± 5.0 + 

- 4.6 

- 0.2 ± 6.5 + 

- 4.5 

- 1.0 ± 8.4 + ^-2 
- 10.7 

352.5 

+ 77 

4.7 ± 8.7 ^ 

- 0.0 

- 3.1 ± 5.0 + 

- 0.8 

1.0 ± 5.6 ++® 
- 0.7 

- 1.7 ± 7.3 + 2 -° 
- 5.9 

- 3.8 ± 12.5 +®-^ 
- 0.0 


TABLE XIV. Cross-section differences for opposite beam helicities in pb-GeV“^ with their statistical and 
asymmetric point-to-point systematic nncertainties for the KinX2 setting, for each bin in (j> (vertical) and 
—t (horizontal). 




if ) ( deg ) 

xb = 0.336 
= 2.161 GeV ^ 
t = - 0.171 GeV ^ 

Xb = 0.342 

Q 2 = 2.190 GeV ^ 
i = - 0.231 GeV ^ 

Xb = 0.343 
= 2.194 GeV ^ 
t = - 0.278 GeV ^ 

Xb = 0.342 
= 2.191 GeV ^ 
t = - 0.324 GeV ^ 

xq = 0.342 

Q 2 = 2.193 GeV ^ 
t = - 0.371 GeV ^ 

7.5 

7.3 ± 6.5 + 

- 8.0 

14.4 ± 10.1 + 

- 6.4 

- 

- 

- 

22.5 

1.7 ± 6.7 

- 1.8 

±03 

32.4 ± 8.9 

- 3.2 

- 

- 

- 

37.5 

7.6 ± 6.1 + 

- 0.4 

16.2 ± 7.9 ^ 

- 0.0 

16.8 ± 8.7 + 

- 7.4 

5.3 ± 11.9 + 2 -® 
- 3.5 

±56 

- 23.9 ± 23.3 ^ 

- 11.7 

52.5 

22.5 ± 5.6 ^ 

- 0.0 

17.7 ± 6.4 

- 2.7 

11.7 ± 7.4 +®'® 
- 1.6 

16.5 ± 9.0 +°-° 
- 4.3 

7.7 ± 12.1 + 

- 0.0 

67.5 

25.8 ± 5.3 ^ 

- 0.0 

20.2 ± 5.8 ^ 

- 1.6 

12.5 ± 6.5 + 

- 4.6 

17.9 ± 8.0 + 

- 5.4 

5.9 ± 8.6 +°-2 

- 4.5 

82.5 

14.9 ± 5.1 

- 0.4 

8.7 ± 5.0 

- 3.4 

17.9 ± 5.4 +°-^ 
- 0.9 

±45 

13.0 ± 5.5 ^ 

- 0.7 

0.2 ± 6.1 + ^'^ 
- 1.3 

97.5 

14.2 ± 4.5 

- 0.0 

21.1 ± 4.8 

- 0.8 

13.5 ± 4.6 + 2-4 
- 0.0 

12.5 ± 4.9 +®-° 
- 0.2 

12.7 ± 5.3 + 

- 3.4 

112.5 

23.3 ± 4.2 ^ 

- 2.6 

±14 

18.6 ± 4.0 ^ 

- 0.6 

11.3 ± 4.2 +°-2 
- 2.0 

1.5 ± 4.2 +^-® 

- 0.0 

- 3.2 ± 4.3 + ^'^ 

- 1.0 

127.5 

11.3 ± 4.2 ^ 

- 3.8 

11.5 ± 3.6 ^ 

- 4.6 

2.3 ± 3.3 + ^'2 
- 0.8 

9.5 ± 3.9 +°-^ 
- 0.7 

5.1 ± 3.6 +°-^ 

- 2.0 

142.5 

7.0 ± 3.9 ^ 

- 0.2 

10.7 ± 3.5 ^ 

- 0.7 

±17 

5.7 ± 3.4 ^ 

- 0.0 

7.8 ± 3.7 +°-° 
- 1.1 

4.5 ± 3.5 +°-® 

- 1.9 

157.5 

- 3.1 ± 3.7 + 

- 0.0 

6.4 ± 3.5 

- 0.9 

4.1 ± 3.5 + °-2 
- 3.0 

7.2 ± 3.6 +°-° 
- 1.0 

3.8 ± 3.7 + 

- 0.4 

172.5 

4.5 ± 3.7 

- 0.3 

1.4 ± 3.4 + 

- 0.0 

0.8 ± 3.6 + 

- 0.9 

0.8 ± 3.4 + 

- 1.9 

- 0.2 ± 3.9 + 

- 0.0 

187.5 

- 1.7 ± 3.5 

- 2.5 

- 1.2 ± 3.4 + 2.2 

- 0.5 

2.5 ± 3.6 + 

- 1.9 

- 0.2 ± 3.4 + 

- 1.0 

8.5 ± 3.6 +^-^ 

- 0.3 

202.5 

- 2.1 ± 3.6 ^ 

- 1.3 

1.4 ± 3.6 + 

- 1.1 

- 3.3 ± 3.3 + 

- 0.1 

- 9.4 ± 3.8 + 2 -® 
- 0.0 

0.5 ± 3.8 + 

- 3.6 

217.5 

- 9.2 ± 3.9 ^ 

- 0.7 

- 9.5 ± 3.7 +°-° 
- 2.2 

±14 

- 6.3 ± 3.6 ^ 

- 0.0 

- 7.3 ± 3.4 + +2 
- 0.0 

- 4.1 ± 3.4 + 2 '® 
- 0.2 

232.5 

- 13.5 ± 4.1 

- 0.0 

- 3.8 ± 3.8 + 

- 2.1 

- 3.4 ± 3.7 + 

- 0.5 

- 3.9 ± 3.8 + 

- 0.8 

±13 

- 7.9 ± 3.8 ^ 

- 0.2 

247.5 

- 6.8 ± 4.1 + 2 -'^ 
- 0.0 

- 8.4 ± 4.0 +°-° 
- 1.8 

- 3.6 ± 3.8 + °'2 
- 0.9 

- 0.6 ± 3.8 + 

- 0.0 

- 4.9 ± 4.2 + 2 + 

- 0.5 

262.5 

- 22.7 ± 4.5 ^ 

- 0.8 

- 11.2 ± 4.5 +°-® 
- 0.7 

- 19.1 ± 4.3 + 

- 5.1 

- 3.0 ± 4.8 + °'’^ 
- 2.0 

- 0.9 ± 4.4 + 2 -® 
- 0.1 

277.5 

±04 

- 19.9 ± 4.8 ^ 

- 1.4 

- 22.1 ± 4.8 + 2 -^ 
- 0.8 

- 14.7 ± 5.0 + 

- 0.0 

- 19.4 ± 5.6 +°-° 
- 3.1 

- 4.8 ± 5.7 +°-° 
- 2.5 

292.5 

- 17.6 ± 5.4 + 2 -^ 
- 0.5 

- 21.5 ± 5.4 + 

- 1.6 

- 16.1 ± 5.5 +°-® 
- 2.7 

4 - 1 ^ 

- 20.0 ± 7.1 ^ 

- 2.7 

- 15.4 ± 7.7 +°-° 
- 3.8 

307.5 

- 23.2 ± 5.5 ^ 

- 4.8 

±79 

- 27.4 ± 6.3 ^ 

- 0.0 

- 3.2 ± 6.9 + 

- 3.4 

- 21.8 ± 7.8 + 2-2 
- 1.8 

- 9.7 ± 10.6 + 2 -^ 

- 1.5 

322.5 

- 19.2 ± 5.9 + 

- 0.0 

- 17.0 ± 7.2 + 

- 2.9 

- 2.6 ± 8.1 + 

- 0.0 

- 26.5 ± 12.3 + 

- 0.0 

- 19.9 ± 17.6 + 

- 25.8 

337.5 

- 4.7 ± 6.2 

- 0.9 

5.9 ± 8.3 +°-° 
- 5.4 

- 

- 

- 

352.5 

- 4.1 ± 6.5 + 

- 3.4 

6.6 ± 9.4 + 2 -® 
- 2.0 

- 

- 

- 


TABLE XV. Cross-section differences for opposite beam helicities in pb-GeV“^ with their statistical and 
asymmetric point-to-point systematic nncertainties for the KinX3 setting, for each bin in (j> (vertical) and 
—t (horizontal). 
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FIG. 24. (Color online) Unpolarized cross sections for Kin2. Each t-bin corresponds to slightly different av¬ 
erage {xb, Q^) values; their range is indicated in the legend, their specific values are listed in the data tables. 
Error bars are statistical only. The light blue area represents the point-to-point systematic uncertainties 
added linearly to the normalization error. The KMlOa model along with its modified version (including the 
TMC effects) are shown as dotted blue and solid green curves, respectively. The Bethe-Heitler contribution 
is represented as a dashed red line. 


[5] A. V. Radyushkin, Phys. Rev. D56, 5524 
(1997), hep-ph/9704207 

[6] A. V. Radyushkin, Phys. Lett. B380, 417 
(1996), arXiv:hep-ph/9604317 [hep-ph] 

[7] J. C. Collins, L. Frankfurt, and M. Strik- 
man, Phys. Rev. D56, 2982 (1997), arXiv:hep- 
ph/9611433 [hep-ph] 

[8] X. Ji and J. Osborne, Phys. Rev. D58, 094018 
(1998), hep-ph/9801260 

[9] A. V. Belitsky, D. Mueller, and A. Kirch- 
ner, Nucl. Phys. B629, 323 (2002), hep- 
ph/0112108 

[10] M. Burkardt, (2007), arXiv:0711.1881 [hep- 
ph] 

]11] O. V. Teryaev, (2005), arXiv:hep-ph/0510031 
[hep-ph] 


[12] I. V. Anikin and O. V. Teryaev, Phys. Rev. 
D76, 056007 (2007), arXiv:0704.2185 [hep-ph] 

[13] I. V. Anikin and O. V. Teryaev, Fizika B17, 
151 (2008), arXiv:0710.4211 [hep-ph] 

[14] M. Diehl and D. Y. Ivanov, (2007), 
arXiv:0712.3533 [hep-ph] 

[15] M. V. Polyakov and C. Weiss, Phys. Rev. D60, 
114017 (1999), arXiv:hep-ph/9902451 [hep-ph] 

[16] A. Airapetian et al. (HERMES), Phys. Rev. 
Lett. 87, 182001 (2001), hep-ex/0106068 

[17] S. Stepanyan et al. (CLAS), Phys. Rev. Lett. 
87, 182002 (2001), hep-ex/0107043 

[18] C. Adloff et al. (HI), Phys. Lett. B 517, 47 
(2001), hep-ex/0107005 

[19] S. Chekanov et al. (ZEUS), Phys. Lett. B573, 
46 (2003), hep-ex/0305028 





































41 



FIG. 25. (Color online) Unpolarized cross sections for Kin3. Error bars are statistical only. The light blue 
area represents the point-to-point systematic uncertainties added linearly to the normalization error. The 
KMlOa model along with its modified version (including the TMC effects) are shown as dotted blue and 
solid green curves, respectively. The Bethe-Heitler contribution is represented as a dashed red line. 
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FIG. 26. (Color online) Cross-section differences for opposite beam helicities for Kinl. Error bars are 
statistical only. The light blue area represents the point-to-point systematic uncertainties added linearly 
to the normalization error. The KMlOa model along with its modified version (including the TMC effects) 
are shown as dotted blue and solid green curves, respectively, except for the first f-bin which is outside the 
prescribed range of this model. 


polarized Protons. 

[32] A. Camsonne, C. E. Hyde, M. Mazouz, 

C. Munoz Camacho, J. Roche, et al. (the Hall 
A DVCS Collaboration), “E08-025,” (2008), 

JLab Experiment E08-025, Measurement of 
the Deeply Virtual Compton Scattering off the 
Neutron. 

[33] F.-X. Girod et al. (CLAS Collaboration), 

Phys. Rev. Lett. 100, 162002 (2008), 

arXiv:0711.4805 [hep-ex] 

[34] S. Chen et al. (CLAS Collaboration), Phys. 
Rev. Lett. 97, 072002 (2006), arXiv:hep- 
ex/0605012 [hep-ex] 

[35] E. Seder, A. Biselli, S. Pisano, and S. Niccolai 
(CLAS Collaboration), Phys. Rev. Lett. 114, 
032001 (2015), arXiv: 1410.6615 [hep-ex] 

[36] S. Pisano et al. (the CLAS Collaboration), 
(2015), arXiv:1501.07052 [hep-ex] 


[37] J. Roche et al., (2006), arXiv:nucl-ex/0609015 
[nucl-ex] 

[38] A. Bacchetta, U. D’Alesio, M. Diehl, and G. A. 
Miller, Phys. Rev. D70, 117504 (2004), hep- 
ph/0410050 

[39] M. Diehl, T. Gousset, B. Pire, and J. P. 
Ralston, Phys. Lett. B411, 193 (1997), hep- 
ph/9706344, 

[40] A. V. Belitsky and D. Mueller, Phys. Rev. 
D82, 074010 (2010), arXiv: 1005.5209 [hep-ph] 

[41] P. Berlin, C. E. Hyde, R. Ransome, F. Sabatie, 

et al. (the Hall A DVCS Collaboration), “EOO- 
110,” (2000), JLab Experiment EOO-110, 

Deeply Virtual Compton Scattering at 6 GeV. 

[42] Formally, Deeply Virtual Compton Scattering 
or DVCS refers only to the sub-process 7 *p —>■ 
7 p. However, DVCS is often used more loosely 
in the literature to name the photon electropro- 













































43 



FIG. 27. (Color online) Cross-section differences for opposite beam helicities for Kin2. Error bars are 
statistical only. The light blue area represents the point-to-point systematic uncertainties added linearly 
to the normalization error. The KMlOa model along with its modified version (including the TMC effects) 
are shown as dotted blue and solid green curves, respectively. 
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FIG. 28. (Color online) Cross-section differences for opposite beam helicities for Kin3. Error bars are 
statistical only. The light blue area represents the point-to-point systematic uncertainties added linearly 
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are shown as dotted blue and solid green curves, respectively. 
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FIG. 29. (Color online) Unpolarized cross sections for KinX2. Error bars are statistical only. The light 
blue area represents the point-to-point systematic uncertainties added linearly to the normalization error. 
The KMlOa model along with its modified version (including the TMC effects) are shown as dotted blue 
and solid green curves, respectively. The Bethe-Heitler contribution is represented as a dashed red line. 
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FIG. 30. (Color online) Unpolarized cross sections for KinX3. Error bars are statistical only. The light 
blue area represents the point-to-point systematic uncertainties added linearly to the normalization error. 
The KMlOa model along with its modified version (including the TMC effects) are shown as dotted blue 
and solid green curves, respectively. The Bethe-Heitler contribution is represented as a dashed red line. 
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FIG. 31. (Color online) Cross-section differences for opposite beam helicities for KinX2. Error bars are 
statistical only. The light blue area represents the point-to-point systematic uncertainties added linearly 
to the normalization error. The KMlOa model along with its modified version (including the TMC effects) 
are shown as dotted blue and solid green curves, respectively. 
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FIG. 32. (Color online) Cross-section differences for opposite beam helicities for KinX3. Error bars are 
statistical only. The light blue area represents the point-to-point systematic uncertainties added linearly 
to the normalization error. The KMlOa model along with its modified version (including the TMC effects) 
are shown as dotted blue and solid green curves, respectively. 
































